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Abstract. In this work we obtain strong unique continuation results for 
variable coefficient second order parabolic equations. The coefficients in 
the principal part are assumed to satisfy a Lipschitz condition in x and a 
Holder C's condition in time. The coefficients in the lower order terms, i.e. 
the potential and the gradient potential, are allowed to be unbounded and 
required only to satisfy mixed norm bounds in scale invariant L^L^ spaces. 



1. Introduction 

The evolution of the understanding of the strong unique continuation prob- 
lem for second order parabolic equations mirrors and is closely related to the 
corresponding strong unique continuation problem for second order elliptic 
equations. Consequently, we begin with a brief overview of the latter problem. 

To a second order elliptic operator = dig^Wj and potentials V, W in 
we associate the elliptic equation 

(1.1) - AyU = WVu + Vu 

Given a function u G Lf^^{W^) and xq G we say that u vanishes of infinite 
order at Xq if there exists R so that for each integer N we have 

(1.2) / lul"^ dx < c^r^^ , r<R 

JB{xo,r) 

The elliptic strong unique continuation property ESUCP has the form 

Letu he a solution to (11.11) which vanishes of infinite order ESUCP 
at Xq. Then u{x) = for x in a neighborhood of Xq. 

ESUCP type results go back to the pioneering work of Carleman [5] in 
dimension n = 2, later extended to higher dimension in by Aronszajn and 
collaborators [3], [Ij. Their results apply to Lipschitz metrics g but only mildly 
unbounded potentials V and W. A key ingredient in their approach was 
to obtain a class of weighted estimates which were later called Carleman 
estimates. The simplest Carleman estimate has the form 

|||x|~^n||i2 < II |a;p^^AM||i2 
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and holds uniformly for r away from — h A^). This restriction is related 
to the spectrum of the spherical Laplacian. 

Adding some extra convexity to the \x\~'^ weight makes the above estimate 
more robust and allows one to also use it in the variable coefficient case. The 
role played by the convexity was further clarified and explained by Hormander 
[12] , |13j . who introduced the pseudoconvexity condition for weights as an 
almost necessary and sufficient condition in order for the Carleman estimates 
to hold. 

The problem becomes more difficult if one seeks to work with unbounded 
potentials V in or near the scale invariant L2 space. There the Carle- 
man estimates are insufficient. Instead the key breakthrough was achieved 
in Jerison-Kenig [15], where the Carleman estimates are replaced by 
estimates of the form 

IIIXI^'^UII 2n < II Io^I^'^AmII 2n 

II I I "ltt^ ~ II I I iil'^ 
Relevant to the present paper is also the alternative proof of this result which 
was given by Jerison [13], taking advantage of Sogge's [25] spectral projection 
bounds for the spherical Laplacian. In the case of operators with smooth 
variable coefficients Carleman estimates were first obtained by Sogge 



Working with gradient potentials in the scale invariant space W E intro- 
duces an added layer of difficulty. There not even the Carleman estimates 
can hold. Wolff's solution to this in |29| is a weight osculation argument, which 
allows one to taylor the weight in the Carleman estimate to the solution u, 
producing estimates of the form 

||e-^'^(^)u|| 2n + ||e-"'^(^)VrVM|| 2n < lle^'^^^^Awll 2n 

where the choice of depends on both u, W and r. 

Finally, the authors's article [17] combines the ideas above into a nearly 
optimal scale invariant ESUCP result for the elliptic problem, with (i) a 
Lipschitz metric g, (ii) an potential V, and (iii) an almost gradient 
potential W. The present paper is the counterpart of [IT] for the parabohc 
strong unique continuation problem. 

We consider the second order backwards parabolic operator 

(1.3) P = dt + dkg'\t,x)di 

in M X and potentials V,Wi,W2- To these we associate the parabolic 
equation 

(1.4) Pu = Vu + WiV^u + V^iW2u) 

Given a function u G Lf^^ and (^052:0) G M x we say that u vanishes of 
infinite order at (to, Xq) if there exists R so that for each integer N we have 

(1.5) [ I \u\^ dxdt < cl^r'^^ , r<R 

Jo JB{xo,r) 



Alternatively we may only require that x —>■ u{to,x) vanishes of infinite order 
at xq, i.e. 

(1.6) / \u{to,x)\'^dx <c%r^^, r<R 

The two conditions f 1 1.5 1) and fll.6p are largely equivalent provided that the 
coefficients g''^ have some uniform regularity as t — >■ to. However, our as- 
sumptions in this article are not strong enough to guarantee this, therefore we 
consider the two separate cases. 

Now we can define the strong unique continuation property SUCP(I) : 

Let u be a solution to ljil-4\ ) which vanishes of infinite order SUCPfl) 
at {tQ,XQ). Then u{tQ,x) = for x in a neighborhood of Xq. 

and the slightly stronger variant 

Let u be a solution to ^1-4^ so that x^u{tQ, x) vanishes of infinite jg-yCPfll) 
order at xq. Then u{to,x) = for x in a neighborhood of xq. 

The study of unique continuation for parabolic equations began with early 
work of Mizohata [2T] and Yamabe [301, followed by Saut-Scheurer ^3]; 
Carleman estimates were first obtained by Sogge 

The study of the parabolic strong unique continuation problem began with 
work of Lin [20j who considered SUCP(II) for the heat equation with W = 
and V bounded and time independent. This continued with work of Chen |6] 
and Poon [22]. Fernandez [11], and Escauriaza, Fernandez and Vessella [7] con- 
sidered SUCP(II) under various assumptions on the coefficients and pointwise 
bounds for W = and V. It is a consequence of Alessandrini and Vessella [1] 
that SUCP(I) andSUCP(II) are equivalent under weak assumptions on the 
coefficients, and they derived SUCP(II) in [2] for bounded W and V. 

The article of Poon [22] contributed to clarifying the correct form of the 
Carleman estimates for the parabolic strong unique continuation problem 
in Escauriaza and Fernandez's work p]. In the simplest form, these have the 
form 

||t~^~2e~ 8^M||i2 < \\t~"''^2e~'st(^dt + A)u\\l2 
and hold uniformly with respect to r away from {2n + N)/4. This restriction 
is connected with the spectral properties of the Hermite operator. 

The spectral projection bounds for the Hermite operator were indepen- 
dently obtained by Thangavelu [28j and Kharazdhov [16]; see also the sim- 
plified proof in the authors's paper [TH]. These bounds were essential in the 
proof of Carleman inequalities for the heat operator of Escauriaza [5] and 
Escauriaza and Vega [10] which yield SUCP(I) when g = In, W = and 
V e L^L°^ + L^L""/^. 

Our aim is to prove that SUCP(I) respectively, SUCP(II) hold under 
sharp scale invariant assumptions on the metric g and conditions on the 
potentials V and Wi, W2- The contribution of this work is comparable to [T7] 
for the elliptic problem: We study almost optimal conditions on 
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(1) the coefficients g 

(2) the potential V 

(3) the gradient potentials Wj 

The combination of rough variable coefficients and conditions on the po- 
tential seems to be new. Also, to the best of our knowledge this is the ffist 
result on unique continuation for parabolic problems under conditions on 
the coefficients of the gradient term. 

For simplicity we always assume that to = 0, Xq = 0. For SUCP(I) it is 
naturaQ to consider a larger class of operators P which have the form 

(1.7) P=^ + dkg'^di + ^d'^di + dd'^^ + ^e''^ 
^ ^ dt ^ t t t t 

where {g^''), {d^^) and (e^') are real valued and {g^^) and (e*^') are symmetric. 
Then simple scale invariant assumptions for the coefficients would be 

(1.8) iMlUoo + ||(t + x'Y^dAU^ + II WIU- < 1 

Here and in the sequel d stands for a generic coefficient of the form g^^ — 
d^^ and e^K For V and Wi^2 we could consider conditions of the form 

(1.9) ||V||^l^oo_^iooJ^n/2 < 1, 

(1.10) ||VFl,2i|L2ioo+iooin < 1. 

Here and in the sequel we use the notation L^L^ = L^L^. 

The situation is however more complex and we may take (11.81) to (11.101) 
only as guidelines. We will have to strengthen (11.81) to include some dyadic 
summability. on the other hand we are able to weaken the time differentiability 
to a Ca Holder condition on small time scales. 

We are also able to slightly weaken (II. 9p almost to uniform bounds on 
dyadic sets. However, we are unable to use mixed norms for Wi and and 
we restrict ourselves to a summable norm in dyadic sets. 

To state our assumptions on g,V, Wi and W2 we consider a double infinite 
dyadic partition of the space, 

00 00 

(1.11) M+ X = y [jAij 

i=—oo j=0 

where 

(1.12) A,j = {(t,x) e M+ X I e-^'-^ < t < e-^\ < 1 + 2|x|r^ < e-^'+^j. 
Consider the subset of indices 

(1.13) A={{t,j):j <2z + 2} 
defining a partition of the cylinder 

Q= [0,1] X 5(0,1) 



""^This becomes clearer later on after a change of coordinates and conjugation with respect 
to a Gaussian weight 
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Define 

(1.14) A{t) = eA:4i> lnr + l,j < ^lnr + 2} 

which corresponds to a partition of the cut parabola 

Qr = {it,x) : |a;|^ < rt < 1}. 




Figure 1. The cut parabola 

We also consider a decomposition of Q into dyadic time slices 

A=[e-^'-\e-^']xB{0,l) 
and a similar partition of the cut parabola Qr into the sets 

Al = Ain Qr 

Given a function space X and 1 < g < cx) we introduce the Banach spaces 
["^{AjX) with norms 



WW" = V 



q 

X{A,,)- 



In a similar manner we define the spaces 1°°{A,X). 

Within the sets Aij we define the modulus of continuity [rriij) in time 

mij{p) = e^'p + e5(2*"-^')pi 

and denote by C'™'-' the space of continuous functions with finite seminorm 

,, ,, \u{ti,x) - U{t2,x)\ 
■U ^ij = sup 7—] ^- 

^* t,M,. m*^(|ti-t2|) 
For the reader's convenience we note that within A^j we have e^* ~ and 

For the coefficients of the operator P in fll.7p we change the condition (11.81) 

to 

(1.15) sup IMIU-(A,,) + e^-2lM||Lip^(A,,^) + < 1. 

^ Mr) 

where we note that e-'"^* ^ (t + x^)^ in Ajj. 
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The pointwise bound for g — In in fll.lSp . namely 

(1-16) sup ll^f - In\\lHA{T),L^) < 1 

T 

is not really needed for our results. It can be always obtained from the other 

bounds after a change of coordinates. This is discussed in the appendix. 

1 

The assertion ( ll.lSp is satisfied for g G Lip^. flC/ provided that 5^(0, 0) = In- 

1 

Indeed by scaling we may assume that the Lip^. fl norm is small therefore 
it suffices to compute 

i>lnT j<ln{T)/2+2 

For the potential V we consider: 

(1.17) IIV^II 

\\V\\i--(A,L^L^+L^L^) '^'^ forra = 1 
where p' in the second line is the dual exponent. In addition we require that 

(1.18) sup II Xi'^^ II LiL^+L-L^/a < 1 > 2 

with the obvious modifications for n = 1,2, where Xi is the characteristic 
function of the set 

{{t,x) : e-*'-^ < t < e-^\ t-'^^\x\ < ^}. 

Both (fTTTD and ffLTSj) are fulfilled if V" G L^L"^ + with small norm. 

Finally for the gradient potentials Wi_2 we introduce the summability con- 
dition with respect to time slices 

(1.19) sup E l|W^l,2|U"+2(A[) < 1. 

T 

t 

As a consequence of this we note the uniform bound 

(1.20) SUp||W^i,2|Un+2(A,) < 1. 

i 

Now we can state our main results. 

Theorem 1. Let P be as in (11.71) with coefficients satisfying fll.lSp . As- 
sume that the potentials V and Wi^2 satisfy (I1.17p . (ll.lSp and (ll.lQp . Then 
SUCP(I) holds at (0,0) for functions u satisfying (11. 4p . 

It is part of the conclusion that the trace of u at t = exists near x = 0. 
The assumptions on the operator seem to be too weak to imply existence of a 
trace in general. More precisely shall prove 

||M(t, 011^2(5(0,1/8)) < 

for some 5 > 0. 
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The C5 Holder regularity in time for the metric g seems so be new, improv- 
ing the C2 Holder regularity in [9]. It is not clear to the authors whether this 
condition is optimal or not. 

We may replace the assumptions by stronger translation invariant assump- 
tions, 

(1-21) \\9\\Lip. + \\9\\l<l 

(1.22) \\V\\lil°-+l°-l^/2 < 1 

(1-23) E 11^1.2 

i 

Then we also obtain a stronger conclusion. 

Theorem 2. Let P he as in (11.31) with coefficients as in (ll.2ip . Assume that 
the potentials V and W satisfy fll.22p respectively (11.231) . Then SUCP(II) 
holds at (0,0) for functions u satisfying (II. 4p . 

We remark that the condition (ll.2ip is really too strong, and that with some 
additional work (see Remark l2.3l) one can bring it almost to the level of (ll.lSp . 
Precisely, it suffices to replace (11.151) by 

(1.24) sup IMIU-(A.,) + e^-'iMI|Lip,(A.,) + \\d\\ „2 « 1. 

where the slightly stronger time continuity modulus mf^ is given by 

m2.(p) = e^^-2V + e5(2^-V 

However, we cannot keep the additional terms in (II. 7p . because we need to 
be able to meaningfully talk about the trace of the solution at time t = 0. 

Both theorems are consequences of quantitative estimates, which also imply 
weak unique continuation under the assumptions of Theorem [2) 

Let u be a solution to for which n(to, •) vanishes in 

the closure of an open set. Then {uifo, .) vanishes in a UCP 

neighborhood of the closure. 

If u satisfies the assumptions and vanishes in an open set U, then it vanishes 
in the time slices t = to in an open neighborhood of the closure of UtQ = {{x : 

{t,x)eU}. 

Theorems [1] and [2] are nontrivial consequences of a Carleman inequality. To 
state a first version of the Carleman inequality we introduce an additional fam- 
ily B{t) of sets which is a partition of the cylinder [0, r~^) x B{0, 1), consisting 
of 

(1.25) Aij, Inr < 4z < r^/^0 < j < lnr/2 + 2, 

(1.26) [e-^'-\ e-^'] x 5(0, e-^W^/^), 4^ > 
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(1.27) Aij, Inr < Ai, lnr/2 <j< 2i. 

This partition is coarser than the partition of the same cyhnder into the sets 
Aij. This is the reason why we need the assumption fll.lSp . More precisely 
Assumptions (11.171) and (11.181) imply 



(1.28) Ir II/°°(b(t),lil°°+l°°l"/2) 1- 

Theorem 3. Let tq ^ 1, £ > and P as in ( II. 7p with coefficients satisfying 
( ll.lSp . Suppose that Wi^2 satisfY (I1.19P with constants depending on e and 
tq. Then there exists C > such that for all r > tq the following is true: 
Suppose that v G L'^{H^) is compactly supported in [0, 8r-i) X 5(0,2) and 
that it vanishes of infinite order near (0,0). Then we can find a function 
G C~([0,8r-i] X B{0, 1)\{0,0}) and h G (:7°°(M+) which satisfy 

(1.29) T<h'<{l+e)T 



(1.30) 
such that 



{x,t)- [h{-\nt)- 



x"^ 



for n > 3, respectively 

\\<^^^\\lHB{T),L°-L^nLP'Li') < C\\e'^{P + WiV + VW2)v\\i2^i3(^^)^LiL2+LPLi), 

forn = 2, i + i = |, 2 < and 
for n = 1. 

The statement of the Carleman inequahty is involved for several reasons. 
The weight t~'^e~^^^ which works for the constant coefficient case, has to be 
modified so that it has more convexity in order to handle variable coefficients, 
spatial localization and the gradient potential. However the polynomial growth 
in time (imposed by the assumption of vanishing of infinite order) limits the 
available amount of convexity; this is the origin of the summability in (I1.16p . 
( I1.20p . and to a lesser extend of (I1.18p . 

If ly = then the Carleman inequality holds for a large explicit class 
of weights e"^. This cannot be possibly true for general gradient potentials. 
Instead, we are only able to prove that there exists some weight function 0, 
which now depends on r, u and W, for which the uniform Carleman inequality 
holds. This strategy goes back to the seminal work of T. Wolff f29j and has 
been used by the authors for the elliptic problem [T7] . 

The partition A^j is much finer than the dyadic decomposition in t only, 
which would correspond to the dyadic decomposition in the elliptic case. We 
are only able to localize the estimates to the sets Aij if we make the weight 
function sufficiently convex. We can do this for many A^j, but not for all of 



them. The sets fll.26p correspond directly to the assumption (11.181) . We need 
to have control of the L^L°° + L°°L^/'^ norm of V in sets which are not smaller 
than those of the partition in fOSD . fOGj) and ffr27l) . 

We have stated Theorem [3] in a simpler form which suffices to derive Theo- 
rem [T] and Theorem [2l However the full estimate we prove is stronger in that 
it also contains precise bounds. These are essential for the localization and 
perturbations techniques we use. 

The strategy of the proof is the same as in [T7]: 

(1) We construct families of pseudoconvex weights and derive Carle- 
man inequalities. The convexity of weights determines the space-time 
localization scales and the admissible size of perturbations. 

(2) We enhance the above Carleman inequalities to include esti- 
mates. Due to the localization it suffices to do this in small sets. 
This allows us to use perturbation arguments starting from the case of 
the heat equation with the weight t~^e~'^' 

(3) estimates for the spectral projections to spherical harmonics im- 
ply the Carleman inequalities in the elliptic case. Here spectral 
projection for the Hermite operator play a similar role. 

(4) Finally we include Wolff's osculating argument into the scheme in order 
to handle the gradient potentials. The efficiency of this part depends 
on the flexibility in the choice of the weight functions. 

The complexity of the weights and the Carleman estimates comes mainly 
from the geometry of the classical harmonic oscillator. Orbits are contained in 
a sphere in M^". The projection down in the x space is a ball, where frequency 
variables have a different behavior in radial and angular directions and near 
the boundary of the the ball. It turns out that our analytic estimates reflect 
these features. 

2. Proof of Theorem □ and [2] 

In this section we prove Theorem [T] and [2] assuming Theorem [31 The relation 
between Carleman estimates and unique continuation is fairly straightforward 
in the elliptic case. In the parabolic situation the argument is less direct due 
to the more complex geometry of the level sets of the weight functions. 

It is a standard consequence of a localized energy inequality that for the 
parabolic equation (11. 4p u{t) and its gradient can be controlled by norms 
of u. 

Proposition 2.1. Let n > 3 and suppose that v solve the parabolic equation 

vt + dka^^div = WiVv + V{W2v) + Vv 

on the space-time cylinder Q = [0, 2] x B{0, 2) with a G Lip uniformly elliptic 
and 

||W^1,2||l"+2{Q) + ||^||l1L°°+L°°L"/2 ^ 1 
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t 

2/r 
1/r 

to 

l/16r 
l/8r 



Es 26 1/4 12^ 
Figure 2. The sets E^, F^^^* and F;'^* 

r/ien 

sup \\v{t)\\L2(^B{0,l)) + \\'^xV\\l2{[0,1]xB{0,1)) < 11"^^ II L2(q) . 
0<t<l 

If n = 2 then the same statement is true with L°°L^/'^ replaced by L^L'^ with 
I < p < oo, 1 < g < oo and ^ + ^ = 1. Similarly, if n = 1 we have to replace 
It by L^L°^. 

Given our assumptions fll.lSp . fll.l7p . (11. 181) . fll.l9p and fll.20p we can apply 

Proposition 12.11 rescaled in sets of the form [to, 2to] x -B(xo,V )' which are 
subsets of the Aij. Summing up with respect to such sets contained in a 
parabohc cube 

Qr = [0,r']xB{0,r) 
we obtain the following consequence. 

Corollary 2.2. The following estimate holds under the assumptions of Theo- 
remUl 

P^m||l-L2(Q,) + \\t^Vu\\L2(^Q,) < ||m||l2{Q2,) 

Proof of TheoremUl We choose r > tq and < 5 <^ r^^/^ and introduce the 
sets 

Es =([0,25'] X B{0,25))\{[0,5^] x 5(0,5)) 
=([0,2/r] X 5(0,2))\([0,l/r] x 5(0,1)) 
F™* =[l/32r, l/16r] x 5(0, 1/4) 

Our strategy will be to truncate u in Es and F^^* and to apply Theorem [3] 
to the truncated function in order to obtain a good bound on u in 

Let 77 be a cutoff function supported in [0,2) x 5(0, 2) and identically 1 in 
[0, 1] X 5(0, 1). For S < r-^/' „e define 

vs{t, x) = (1 — r]{t/5'^, x/(5))?7(rt/8, x)u{x, t) 
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which satisfies 

(P + WV)v8 = Vvs -[P + WV, v{t/S^ x/6)]u + [P + WV, r/(rt/8, x)]u. 

The second term on the right hand side is supported in Es and the third one 
in F^^*. 

We apply Theorem [3] to vs- One should keep in mind that the corresponding 
weight (j) depends on 6 but that the bounds we prove are uniform with respect 
to 6. We normalize the function h by h{0) = 0. We have to control the size of 
(j) in the sets Es, F^* and F^'^*. Due to (Ol) we have 

Ts < h{s) < 2ts, s>0 

By (11.301) we obtain a rough polynomial bound in 6 



.|2 



(2.1) < t-2-e-^+^(^+^) < t'/^c{T)S-'^-' in Es. 
Let 

M = sup{e'^(-"^W)e-^+^("+'r) ; (t,x) G F;^*} 

By (11.291) the supremum is attained at a point {to,Xo) with | < StIq < 1 and 
|xo| = 1. A simple computation also shows that 

sup{t-i/2e'^(-i'^(*))e-^+^(^+^) : (t, x) G F;^*} < r'^^M. 
Then M dominates e'^ in F^^*: 

(2.2) < M, t-i/^e* < r^^^M in F,"^*. 
Next we need to bound from below in F™* in terms of M, 

(2.3) inf e"^ > e^^M 

/T'zrit 

T 

To see this we compute for (t, x) G F^"* and sufficiently small e: 

0(t,x) -0(to,xo) > /i(-lnt) -/i(-lnto) + ^-^-2£(r + ^) 

> (--20e)T>-T 

and use dH]) and (Q. 

Theorem [3] apphed to f ^ yields 

(2.4) 

" "P{B(r),L2L^I=^nL°°L2) ~" ''"i2{B(T),L2L^+LiL2) 



2n 



+ ||e'^[P + iyV,77(rt,x)ln|| 
By Holder's inequality we have 

le'^^V^f^ll ^, ll^ll/oorRl'-rWl7'oo4-7'oor"/2-l lle'^fAll _2n 
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Due to the smallness in (11.281) we can absorb this term on the left hand side 
of the inequahty. 

We calculate the first commutator 

f, = [P + WV,7]{t/5\x/5)]u 

= + dug'^'di + + WV)r^{t/5\ x/5) 

+ 26-'g''{dkv){t/6',x/6)diu 
By (EH]) we have 

For the W term we use (11.191) and Holder's inequality. For term involving 
d^'- we bound the L^L"^ norm in terms of an L°°L'^ norm, using (ll.lSp which 
implies that the pointwise bound for rf^' is summable with respect to dyadic 
time regions. For the remaining terms we simply bound the L^L^ norm in 
terms of the norm. This yields 



+ Ud,ri)it/5',x/S)t-^u\ 



L°°L'^i 



Then we can apply a straightforward modification of Corollary 12.21 on the Es 
scale to finally obtain 

||e^[P + ^yV,r/(^/5^x/5)]u|l < ciT)5~^^-^\\u\\LHE,)■ 

Similarly we can estimate the second commutator 

\\e^\P + WV,r]{Tt,x)]u\\ ^ ^ < MT^/'^\\u\\L^F^-t)■ 

Hence by inequality (12.41) we get 

V ' II ''llp(B(T),L2L7J^nL°°L2) ~ II 11^ l^r i ' V / II 11^ Kt^s) 

Within F^"* we have = u. Then by (12.31) we obtain 

(2.6) ||m||l°°l2{f^"*) < r^/^e"^^||M||L2(F^e.t) + c{t)5''^''''^\\u\\l2[Ei,)- 

Also by the vanishing of infinite order the second term tends to zero as 5 ^ 0. 
Hence as 5 — we obtain 

(2.7) \\u\\L°°L-^{Fi^t) < r^/^e"2^||^x||i2(^ext) 
For < t ^ 1 we choose r = -|- to obtain 

||M(i,-)IU2(iJ(0,l/4)) <t~^/V^. 

This completes the proof of Theorem [H □ 
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Proof of Theorem [3 We extend the potentials V and W by zero to negative 
time, and g^\t,x) = g^\0,x) for t < 0. By definition, possibly after rescaling, 
we have u{0, .) G L^{B{0, 2)). We now solve the mixed problem 

(2.8) ut + dkg'^^diu = for t < and |a;| < 2 
with the boundary condition 

u{t, x) =0 if |x| = 2 and t <0 

and the obvious initial condition to obtain an extension of u to negative. The 
heat kernel for (12.81) satisfies Gaussian estimates. In particular we obtain from 
(II. 6p for all positive integers N with a constant possibly differing from (11.61) 

(2.9) r [ \u\^dxdt < clr^"" . 

We seek to prove that the bound (12. 7p still holds in this context. The 
difficulty is that we only know that u vanishes of infinite order at (0, 0) for 
negative time. To account for this we shift the time up, t ^ t + 26. Arguing as 
in the previous proof we obtain (12. 6p with u replaced by u(t + 26, x). Letting 
5 — s> by (12. 9p we obtain (12. 7p and conclude as above. 

Remark 2.3. If one wants to prove Theorem\^ under the weaker assumptions 
on g in (ll.24p then the origin needs to be avoided in the above argument. Hence 
the time translation needs to be accompanied by a spatial translation, namely 

us(t, x) = u(t + 25^, X — 8r5ei) 

This translation places the image of the origin, or better of the cube [0, tS"^] x 
B{0,4:tS), within the region {rt < x^}. But in this region the conjugated 
operator P^, introduced later, is elliptic so only pointwise bounds for g are 
needed for the Carleman estimates. 

□ 

3. BOUNDS IN THE FLAT CASE AND THE HERMITE OPERATOR 

In this section we prove the simplest possible Carleman estimate for the 
constant coefficient backward parabolic equation 

dtu + A^.M = / 

This serves as a good pretext to introduce the class of weight functions which 
is later modified for the variable coefficient case. 

We also describe the change of coordinates which turns the backward par- 
abolic operator into a forward parabolic equation for the Hermite operator 
H . In this way we are able to relate the Carleman estimates for the heat 
operator to spectral information for H. 
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Proposition 3.1. Let u G Li^ with compact support away from 0. Then 

(3.1) \\t-^-'^e-'&u\\L2 < ||r^+5e-|^(9t + A)M||i2 

uniformly with respect to r away from {2n + N)/4. 

Proof. In ]R+ x R" we introduce new coordinates (s, ?/) G M x M" defined by 



(3.2) 
Then 



-4s 



t = e 

1 —2s 

X = y 



— = -Ae'^'— - e~%— — = -e'^' — 
ds dt dx' dy 2 dx 



Atidt + A,.) = -— - 2y— + Ay 



4t^+"/^e-w(5i + A,)t-"/^ew = _|_ + - =; -d, - H =: -Pq 



Hence in the new coordinates our operator becomes 

--2y- 
ds dy 

//I y'^ 
If we conjugate it by t"^' e~~st = e~"'^e~~ we obtain 

^ d 

ds 

where H is the Hermite operator 

H = -Ay + y^ 

Then it is natural to define the new functions 

v{s,y) = e-"^e-^«(e^^ ^e'^y), g{s,y) = e^— ^^^e'^ /(e^^ ^e'^y) 
which are related by 

Pov = g^ {dt + A)u = f 
In the new coordinates, the bound (13. ip becomes 

(3.3) We'^'vh^ < lle^^^Po^lU- 
Denoting w = e^'^^v, we conjugate 

e^^'PoV = e^^'PoC-^^'w = {-ds -H + At)w 
and the above bound becomes 

(3.4) \\w\\l2 < \\{-d,-H + AT)w\\L2. 
Since ds and H — At commute we expand 

\\{-ds - H + Ar)w\\l2 = \\dsw\\l2 + \\{H - 4r)w\\l2 > d{4T,n + Nf\\w\\l2. 

Note the spectral gap, which is essential in order to obtain strong unique 
continuation results. □ 

For later use we also record the following slight generalization of the above 
result. For expediency this is stated in the {y, s) coordinates, i.e. in the form 
of an analogue of (13.31) . 
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Proposition 3.2. Let h be an increasing, convex, twice dijferentiable function 
so that 



d{h\ N) + h" > 



1 



Then 
(3.5) 

ll(i + /,")i/V(%|U2 + 



mill < 1 



i + h' 

for all compactly supported v E . 



L2 



Proof. After the substitution w = e'^^^^v the bound (13.51) becomes 
(3.6) 

(l + /l")^/2 



'i + h"y/'w\\i,+ 



min < 1 



Hw 



< 



H + h'is))w\\l,. 



' 1 + h' 

and we obtain the estimate through expanding the term on the right hand 
side with respect to its selfadjoint and skewadjoint part: 

H + h'{s))w\\l2 =\\dsw\\l, + \\{-H + h')w\\l, + \\{h" 
> [ id{h',Nf + h")\\wis)\\l,ds 



To complete the proof we observe that for each s we have 

\\Hw{s)\\l2 < \\{~H + h')wis)\\L2 + h'{s)\\wis)\\L2 



□ 



4. Resolvent bounds for the Hermite operator 

As seen in the previous section, the spectral properties of the Hermite op- 
erator play an essential role even in the simplest Carleman estimates for 
the heat equation. In this section we take a look at and bounds for its 
spectral projectors and its resolvent. 

The spectrum of if is n + 2N, and its eigenfunctions are the Hermite func- 
tions defined by 

_^ 

Ua = Caidy - y)°'e 2, Hua = (n + 2\a\)ua 

As |a| increases, so does the multiplicity of the eigenvalues. We denote the 
spectral projectors by H^ for A G n + 2N. We consider both the spectral 
projectors and the resolvent of H and obtain both and localized bounds. 

4.1. Weighted bounds. We consider two parameters 

l<d,R<X"^ 

We denote 

Br = {y : \y\ < R}, ,B^^ = {y: \yj\ < d}, j = l,...,n 
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By XR: respectively Xd we denote bump functions in Bji, respectively 5^ which 
are smooth on the corresponding scales. 

Proposition 4.1. The spectral projectors Ux satisfy the localized Li^ hounds 

(4.1) i?-h^||xiinA/|U2 + i?-U-3||x/?vnA/|U2 < 

respectively 

(4.2) rf-^lll^.l^(x^,n,/)|U.<||/IU. 



Proof. The inequality (14. ip is trivial unless i? ^ A 2 . To prove it in dimension 
n = 1 we only need to consider the case when / is a Hermite function, 

f = Iixf = hx 
in which case it follows from the pointwise bound 

X-'^\h'^{x)\ + X^hx{x)\ < Whxh^, \x\ < ^X^ 

In dimension n = 1 (14. 2 p follows by interpolation from (14. ip with R = d. 
This extends trivially to higher dimension by separation of variables. 

It remains to prove (14.10 in higher dimensions. Summing up (14. 2 p with 
d = R over j we obtain the bound 

R-^\\D\hRllxf\\L^<\\f\\h 

For |x| < -R <^ we have \^\'^ ~ A in the characteristic set of H~?R.z, therefore 
the above norm should essentially control the left hand side of (14.10 . For later 
use we prove a slightly more general result, which in particular concludes the 
proof of (1411]) . 

(4.3) X^\v\\l2+X-^\Vv\\l2< \\\D\-^v\\L^ + \\\y\-^v\\L^ + UH-X)v\\^^^^^^^^^i^^ 
Indeed the norm on the right is equivalent to 

\\H^v\\l2 + \\{H + Xy^H - X)v\\l'2 > X^\\v\\l2 + X'^\\H^v\\l2 
In our case we apply (14.30 to t> = R~^XR^xf ■ Then 

lll?/|^^^IU2 < ||nA/||L2 

while 

(H - X)v = 2i?-^V(Vxi?n;,/) + R-'^AxRllxf 

which yields 

\\(H-X)v\\^,.,^^,^,^,2<R-Hf\\L2 

□ 
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To state the corresponding resolvent bounds we define the spaces X2{z) by 

M\x2{z) = i^ + \'^^\)^\Mr^ + UH-z)u\\ i^^ + supd-^\\\Dj\^XdU\\L^ 

^ II j,d 

and the corresponding dual spaces These spaces are larger than the 

corresponding "elliptic" spaces, 

(4-4) ll^lUaC^) ^ kl^ll^lU^ + \\yv\\L^ + \\Vv\\l2 

On the other hand by extending the bound (14.31) to complex A we obtain a 
counterpart of (14. ip . namely 

(4.5) R'^X-^WxrAl^ + R~^>^'Hxr^u\\l^ ^ Il^llx2(^) 
Finally, the result of (14. 2p can be written in the following dual forms 

(4.6) l|nA/||x.(A) < ll/IU^, ||n,/|U. < ||/||x|(A) 

The localized resolvent bounds have the form 

Proposition 4.2. Let n > 2, z ^ C with dist{z,n + 2N) > 1, and 1 < d < 
R<^^z. Then 

(4.7) lkllx.{.)< ll(^-^)/llx5(.) 

where the d component of norms is omitted in dimension n = 1. 

Proof. We first note that the bounds (14. 6 p almost imply (14. 7p up to a loga- 
rithmic divergence. They do imply easily a bound for higher powers of the 
resolvent for z away from the spectrum of H, 

(4.8) UH - ^)-^-Vllx.(.) < (1 + k > 1. 
as well as 

(4.9) \\u\\L2<il + \^z\)-'^\\iH~z)uy.^,y 
Hence it remains to show that 

(4.10) < (1 + \Qz\t^\\uh^ + m- z)uy,^^y 

Using a positive commutator technique we first prove a one dimensional 
estimate. For this we define the one dimensional skewadjoint pseudo differential 
operatoiB 

Qr = ^Op'"{x{yr-')xm-')) 
where x is a mollified signum function which satisfies 

X{x) = ^, \x\<2 

Its properties are summarized in the following 

^As defined the symbol of Q is not smootli at 0. However, any smootli modification in 
the ball {x^ + < r^} will do. 
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Lemma 4.3. a) is bounded in for 1 < p < oo uniformly for r > 1. 

b) Qr is also bounded in X2{z) uniformly with respect to z & C and r > 1. 

c) Qr satisfies the commutator estimate 

(4.11) r-i \D\'^Xru\\h < (1 + l^^uWl, + {{H - z)u, Qu) 

Proof, a) The boundedness is straightforward and is left for the reader. 

b) For the X2{z) boundedness we consider first the d terms, which without 
any loss in generality we can write in the form 

d-^\{d^ + D^)-^Xdu\\L^ 
Since Q is bounded in it suffices to prove the commutator bound 

\\[Qr,{d^ + D^y^Xd]u\\v^ < \\u\\l2 

But this is easily verified using the pdo calculus. 
Next we consider the term 

for which it suffices to prove the commutator bound 

or equivalently, by duality, 

\\[Qr,H]u\\L2 < \\yu\\L2 + ||Vm||l2 + |||z|5m||^2 

This follows again from the pseudodifFerential calculus. 

c) Since Qr is skewadjoint we have the identity 

{{H — z)u, QrU) = {[H, Qr]u, u) + '^z{iQrU, u) 
therefore it suffices to insure that 

(4.12) {[H,Q,]u,u)>\\u\\%^^.^ + 0{\\u\\l.). 

For this we compute the commutator [H, Q] , 

[if, Q] = Op-({e? + yl X{yir~')x{^i\yi\-')}) + 0{1)l2^l. 
= Op'"{2r-'x'{yir~')^iXi^i\yi\-')) + 0(1)^^^^^ 

Now 

r-\xinei + r^)^ < r-'x'iy,r-')^MM-') + 1 
and the conclusion follows by Garding's inequality. □ 

We return to the proof of the proposition. By separation of variables the 
bound (14. lip extends to higher dimensions and gives 



r-'\ 



\Dj\hiu\\h < (1 + P^DII^IIi, + 2^{{H - z)u, Qiu) 



where Ql the higher dimensional analogue of Qr with respect to the j variable. 
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The X2{z) boundedness of Qr also extends easily to higher dimension. Hence 
by Cauchy-Schwartz we obtain 

r-'\\\D,\hiu\\h < (1 + l^^DMh + m- z)u\\x*i.)Mx,i.) 
To conclude the proof of the estimate fl4.10p it remains to show that 

which follows by duality from (14.41) . 

□ 

The final step in the resolvent bounds is to replace the y' derivatives by 
angular derivatives. Let V± = ||| A V be the angular derivative and \D±\2 be 
the corresponding fractional derivative. 

We split the coordinates into y = {yi,y') and use the notation ' for coordi- 
nates and derivatives in the obvious sense. For 1 < d < R < ^/X we define the 
sector 

BR^a = {R<\y,\<2R, \y'\<d} 
and XRd^ bump function in Bji a- Then we define the function space X2{X, R, d) 
by 



+ R-"'x"'\\xRMh + d-"'\\\D^\^: 

and ^2 (A, i?, d) as its dual. 

Lemma 4.4. Suppose that n >2 and 1 < d < R. Then 

\\u\\x2{X,R,d) ~ R'^>^^\XrM\l^ + ^"^A"3||VXi?,d?^||L2 + d-^\\D'\^XR,d^\\L^ 

Proof. Within Bji^d the angular derivatives are close to the y' derivatives, 
namely 

\D^u\ < \D'u\ + 4|Vn|, \D'u\ < \D^u\ + ^\Vu\. 

This implies the corresponding bounds for norms, and the conclusion follows 
by interpolation. □ 

/■ From the above lemma we obtain 



\u\\x2{\,R,d) ~ ll^llx2{A,fi,d) 



Hence, we may replace X2 by X2 in (14.71) and (14. 6p : 
Corollary 4.5. a) For A in the spectrum of H we have 

(4.13) \\U,f\\x,(x,R,d) < II/IIl^, IIHa/IIl^ < \\f\\xux,R,d) 

b ) For z away from the spectrum of H and 1 < d < R < ^z we have 

(4.14) \\{H-z)-'-'f\\x2ix,R,,^ < (1 + I^SzD-'^II/IUka,/?,.), k > 
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4.2. The LP bounds of the resolvent. The L'p bounds for the spectral 
projectors and the resolvent were proved in [TB], (see also [ID])- For the 
sake of completeness we also present them here in a simpler manner following 
the approach in [19] . We refer the reader to the same paper for further results. 
We consider pairs of exponents satisfying 

. ^\ 2 n n 

4.15 - + -=0 

where the range for p is 

(4.16) p > 4 for n = 1, p > 2 for n = 2, p > 2 for n > 3. 

This leads to the following rang^ for q: 
(4.17) 

2n 

q G [2, oo] for n = 1, g G [2, oo) for = 2, g G [2, -] for n > 3. 

The dual exponents are denoted by p' and q' as usual. 

Proposition 4.6. Let q be as in ( 14.171) . Then 

a) The spectral projectors Il\ satisfy 

l|nA|L,'^L2 <1, IIHaIU^^l, < 1, n>2 

(4.18) 

\\^m\li'~,l2 ^ ^ ^ l|nA||L2^L9 < A ^ n = i 

b) For z away from n + N the resolvent {H — z)^^ satisfies 
(4.19) \\(H - z)-^\\^,,^^, < (1 + p^|)i-^, n > 1 



Outline. To revisit the If bounds associated to the spectral projectors we recall 
the approach in |T^. The first step there is to establish pointwise bounds for 
the Schrodinger evolutioE@ 

(4.20) ||e^*^|Ui^z.oo <(sint)-t 

This immediately (see also [I8]) leads to Strichartz estimates for the solution 
to the inhomogeneous equation 

ivt - Hv = g, v{0) = vq 

namely 

(4-21) |K'||LP{[0,27r];L'J) < IkolUa + \\9\\ Lp' i[0,2n];Li') 

where {p,q) are as described in fl4.15p . (14.161) . 

To obtain KT8\i we apply fl4:2TD to v = C'^'^Uxu, which yields L^ 
bounds, and hence by duality and selfadjointness all estimates of ( 14.18^ for 

■^The exponent g = oo is actually allowed in the spectral projection bounds in dimension 
n — 2. However, it is not allowed in any of the resolvent bounds. 

^These bounds are very robust and are in effect established in [T2_ for a much larger class 
of operators 
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n > 2. The case n = 1 can be dealt with directly using the pointwise bounds 
for the Hermite functions. 

We note a consequence of the bounds (14.181) . namely 

(4.22) UH - zr'''\\^,,_^^, < (1 + \Qz\)-^-^-\ n>2, k>l 

which is obtained by interpolating between q = 2 and q = 
Similarly we get 

(4.23) WiH - z)-'\\^,^^, < (1 + \^z\)-V, n>2. 
Then we apply (I4.2ip to 

v{x,t) = x{t)e-'Mx), 9 = x'{t)e-'Mx) + x{t)e~''\H - z)u 
where x is a unit bump function on an interval of size (1 + |532;|)^^. This yields 

Ml. < (1 + \^z\)^\u\\l. + (1 + \^z\)"^-7\\[H - z)u\\^, <\\{H- z)u\\^,. 

concluding the proof of (14.191) for n>2. The case n = 1 is a variation on the 
same theme. □ 

4.3. Combining the estimates. Here we combine the and the U com- 
ponents in the resolvent bounds: 

Proposition 4.7. For z away from n + 2N the resolvent [H — z)^^ satisfies 
(4.24) 

\\{H - z)-l^,'^x,(^?.Ad) < (1 + l^^^l)'"^ n>2, (n, q) ^ (2, cx)) 
with the obvious modification for n = 1. 

Proof of Proposition Taking into account the bounds (I4.19P and (I4.23p . it 
remains to prove the estimate 

\Mx,in.,R,d) + l^^^l)^"^ ll(^ - ^Ml.' 

+ (1 + \^z\)^^'p\\u\\li + (1 + |$52;|)3||u||i2 

But this follows from (14.120 in the same way as for Proposition 14.21 since the 
operator Q is bounded in L^. □ 

5. LP ESTIMATES IN THE FLAT CASE AND PARAMETRIX BOUNDS 

In this section we begin with the mixed norm L^L'^ Carleman estimates 
in the simplest case, i.e. with constant coefficients and a polynomial weight. 
These were proved in [8] except for the endpoint which was obtained later in 

m ■ 

After a conformal change of coordinates and conjugation with respect to the 
exponential weight the Carleman estimates reduce to proving L^L"? estimates 
for a parametrix K for dt — H -\- t. In this article we need a stronger version 
of these bounds, where we add in localized norms. 

In a simplified form, Escauriaza- Vega's result in fiU\ has the form: 
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Theorem 4. [TO] Let p and q be as above. Then 

2 2 

for all u with compact support in R" x [0, oo) vanishing of infinite order at 
(0, 0) uniformly with respect to 4r with a positive distance from integers. 

One can write the estimate in the (s, y) coordinates using the same trans- 
formation as in Section [3l 

(5-1) ||e"t;||L°°{L2)nLP(L'J) < \W{ds + -f^)^^||Li(L2)+Lp'(L9') 

Setting w = e'^'^v this becomes 

(5.2) lkl|L-(L2)nLf(L«) < 11(5. + ^ - ^)w^IIl1(L2)+Lp'(L9') 

Denoting by 11^ the spectral projection onto the A eigenspace of H we obtain 
a parametrix K for [dt — H + r), 

K{dt + H-t) = I 

where the s-translation invariant kernel of K is 

ir(s) = ^nAe^(^"^h,(,_A)<o 
AeN 

Since w decays at ±oo we have 

w = K{ds + H — t)w 
therefore fl5.2l) can be rewritten in the form 

(5-3) ||-f^/||L°=(L2)nLP(L9) < ||/|lLi(L2)+Lf'(L9') 

The main result of this section is an improvement of (15.21) . namely 
Proposition 5.1. Assume that r is away from n + N and that 

l<d<R<T 

Then 

(5.4) ||-f^/||L°=(L2)nLP(L9)nL2X2(r,R,d) < \\f\\L^L2)+LP'(Li')+L2X*{T,R,d) 

Proof. We work in dimension n > 2; some obvious adjustments are needed in 
dimension n = 1, which is slightly easier. We consider four endpoints: 
A: The L^L^ L°°LF' bound follows easily since the projectors 11;^ are 
bounded. 

B: The L^L^ — > L^U bound. Here it suffices to prove 

\\K{.)f\\^Li < II/IIl2 

Splitting / into spectral projections and using (14.181) we obtain 

II^W/IU.<5^e-i(^-^)*i||n,/|U2 

A 
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For |t| > 1 we can use Cauchy-Schwartz to obtain 

\\mf\\L.nW)<^-'^'^\\f\\L^ 

which suffices for all q. For \t\ < 1 we consider the most difficult case p = 2 
and compute 
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wmn 



2 



L2([„i,i],L<,) 



0<ij \|A-t|ss2* / \\fi-T\^2^+j 



1 1 

2 / \ 2 



0<j 0<i \|A-r|»2' / \|M-r|«2^+J 



L2 



<II/Ili2 

C: The L^L^ — ^ -^2(t, -R, c?) bound for K follows in the same way from 

l|nM||x2(T,/?,d) < I|m||l2- 

D: The L^'L^' + L'^X;{T,R,d) L°^L'^ bound for K is equivalent to the 
L^L^ — s> L^L'^ n X2{t, R, d) bound for K* . By reversing time this is seen to be 
the same as the L^L"^ L^L'^ bound for K. 

E: The Lp'L<?' + L^{X*{R,d)) L^U n L^X2{R,d) bound. Using (PHM 
and fl4.13p directly yields 

\\K(s)\\ 2u 2„ < Vel^("-^)l < s-^e""^ 

II ^ ^ " +X* ( R,d)^L'^ nX2( R,d) ^ ~ 

AeN 

Similarly we obtain 

\\K(s)\\ ^ <s-h-'', \\K(s)\\ <s-h-''' 

II ^ ^"L2^L^nX2(i?,d) ~ ' II V ^ll2,7rR+-^2(«''')^L2 ~ 

Interpolation with the estimate gives 
and 

If p > 2 then the Hardy-Littlewood Sobolev inequality implies 

II^*/IIl.l.<||/IL.'z.,'. 

and 

11-^ * f\\L2X2{R,d) ~ WfWLp'Li'^ \\K * fhpLi < \\f\\L2X*{R,d)- 
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With obvious changes the analysis is similar if = 1, 2. 
It remains to prove the — > type bounds, namely 

\\Kf\\L^X2(R,d) ^ \\f\\v^XT,(R4) = 1' 2) 

respectively 

11^/11 2„ < 11/11 2n (n>2) 

For this, following an idea in [lOj, we consider a dyadic frequency decompo- 
sition in time. By the Littlewood-Paley theory it suffices to prove the bound 
for a single dyadic piece at frequency 2\ namely 

(5.5) \\SAD,)K]\\ 2n < 11/11 2n (h > 2) 

and its one and two dimensional counterpart. 

Taking a time Fourier transform we can write (for / G iS(M")) 

S;k{a)f = s(2-V) V ^—UJ = s{2-^a){H -r- ta)-'f 

^ A — T — la 

A 

therefore by the inversion formula 

{SjK){t)f= j e'"'s{2-^a){H -T -iaY^fda 

e'*"^(s(2-V)(if - r - i(j)-')fda 
da^ 

Hence using the resolvent bounds (14.191) and (I4.22p and the first line for |t| < 
and the second line for \t\ > 2~^ we obtain 

\\SjK{t)\\ .... . 2„ < 



and the similar estimate in one and two dimensions. The bound on the right 
is integrable in t, therefore (15. 5p follows. □ 

6. Modified weights and pseudoconvexity 

The main result of this section. Theorem [5] is a considerable improvement 
of Section [3l The weights t~'^ in Section [3l while easy to use, satisfy merely 
a degenerate pseudoconvexity condition, in the sense that the selfadjoint and 
the skewadjoint parts of the operator in (13. 4p commute. This is in contrast to 
strong pseudoconvexity where one obtains better bounds from the positivity 
of the commutator. A perturbation argument easily implies an Carleman 
estimate for variable coefficients as soon as = + 0(t). However, even arbi- 
trarily small perturbations of g from J„ at t = destroy the pseudoconvexity. 

To obtain results for general variable coefficients we need a more robust 
weight with additional convexity. A good way of doing this is by adding 
convexity in t and by using a weight of the form 
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with a convex function h. Then we obtain for the heat operator the strength- 
ened estimates of Proposition I3.2[ The assumption of vanishing of infinite 
order forces us to work with functions h with at most hnear growth at infinity. 
This in turn hmits the convexity of h, and hence the gain from the convexity 
in the bounds. 

These Carleman inequahties with the weight e"'^^''^*^ are more stable with 
respect to perturbations. They can be obtained for coefficients satisfying 

(6.1) \g - /„| + (t + \x\')\dtg\ + {t+ \x\'y/^\Vg\ < {t + \x\y. 

with suitable functions h. It is not difficult to weaken fl6.ll) almost to our condi- 
tions fll.lSp and (11.161) . This venue was pursued by Escauriaza and Fernandez 

In this paper we seek to obtain Carleman inequalities and also to handle 
LP gradient potentials. Both require good spatial and temporal localization, 
which depends on the strength of the estimates. The weights e^*^"'"^*^ seem 
to be insufficient for this purpose. 

Consequently we consider a larger class of weights of the form 

lnt)+0(a;t-i/2,- ln(t)) 

having some additional convexity in y = xt~^/'^. Here we think of (f) essen- 
tially as a function of y with a milder dependence on s = — Int. Obtaining 
pseudoconvexity is not entirely straightforward because the Hamilton flow for 
the Hermite operator H is periodic so no nonconstant function of y can be 
convex along its orbits. We note that the projection of the orbits to the y 
space are ellipses of size 0(-\/r) where r is the energy, centered at 0. Hence 
we can choose to be convex in y for \y\ -C a/t. We compensate the lack of 
convexity of when \y\ ^ ^Jt by the s convexity of h. To elaborate this idea 
we explain the precise setup. 

Let ^1 be small positive constant. We begin with constants {aij}^ (see (11.131) 
and (I1.14P for the notation) which control the regularity of the coefficient^ 
^ki _ ^ki ^ ^ki ^ki p gjygn fiup as in (ll.lSp . 

(6.2) Sia.j = + e^'^'\\d\\upjA,,) + \\d\\^^.j ^^^^y 
The condition (11.150 guarantee that for all r > 1 

(6.3) - < 1- 

We first adjust the ai/s upward so that they vary slowly and do not con- 
centrate in irrelevant regions. This readjustment depends on the choice of the 
parameter r. 

Lemma 6.1. Let aij be a sequence satisfying (16.31) . Then for each r ^ 1 
there exists a double sequence {£ij)A(T) with the following properties: 



'denoted generically by d here and later 
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(1) For each G A{t) 

(2) We have Sij G 1^{A{t)), 

\\^tj\W{A[T)) < 1- 

(3) The sequence Sij is slowly varying, 

lln^iiji -Ine^ajJ < ^(Ki -i2\ + |ji -J2I), («i,Ji), («2,J2) G ^(r). 

(4) The sequence defined by 

£i= ^ e^ij, z>lnr. 

j.{i,j)&A{r) 

has the following properties 

(6.4) I ln£j^ - In^jJ < ^|i2 - ill, e^ij < e^i, £:i[inr/2]+2 ~ 

(5) For each i > Inr i/iere exists an unique < < [lnr/2] +2 w^/i i/ie 
following properties: 

(6.5) ^ Si. 

and 



(6.6) 
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> e"V-^/2 ^j^-(^) < ^- < [lnr/2] + 2. 



We shall see that j{i) is an important threshold. If j > j{i) then we can 
localize our estimates to the corresponding Aij and even to smaller sets. On 
the other hand, we cannot localize to sets smaller than 

(6.7) B,o= U Ai^. 

Proof. To fulfill the conditions ([I])-(11D we simply mollify the atj, 

[lnT/2]+2 

For the last part of (jlj) we redefine 

p.. — p.. I p-|li-([lnr/2]+2)|~. 

This also increases ii by a fixed factor. 

For ([5]) we begin with a preliminary guess for j{i) which we call jo{i) £ I^^- 
We consider three cases. 

lnr/2 if Ei < r"\ 

Ut) = { -ln(eV/2) if r-i<e,<r-V2, 
if < E,. 
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We define 

(6.8) := max{£,,, 2e-^l^'-^°«l£,}, 

which is still slowly varying because Si is slowly varying. 

We define jii) according to (16. 6p . It is uniquely determined since the se- 
quence Sij is slowly varying compared to e~^T2. Since iij is slowly varying we 
must have iij < Eijl. This allows us to conclude that for j close to jo(^)) the 
second term in (16.81) is larger than the first one, 

= 2e-^l^-^°«l£, for |j - jo(OI < 2. 

If jo(i) = then e^q = 2ei > 2t-^I^ and hence jii) =0. If < jo(0 < lnT/2 
then for |j — jo(OI < 2 we have 

therefore jo(^) — 2 < jii) < jo{i). If jo(0 = lnr/2 then for \j — jo{i)\ < 2 we 
have 

£ij < 2e-5li-io»le-io«^-i/2 

and we arrive again at jo(^) — 2 < In all three cases we have \jo{i) —jif) \ < 
2 therefore (16. 5p holds. We observe that Si ^ ii. □ 

The sequence (£ij)^(r) is used to describe the amount of spatial convexity 
needed in the region Aij, which will be reflected in the construction of below. 
The partial sums Si measure the amount of s-convexity needed in [2,^ + 1]. The 
purpose of part ([5]) above is to correlate the two amounts in a region where 
they have the same strength (where j is close to 

Our weights have the form 

(6.9) ^{s,y) = h{s) + (P{s,y) 
Their choice is described in the next two lemmas: 

Lemma 6.2. Let r and (^j) he as in Lemma \6.1\ Then there is a convex 
function h with the following properties: 

(1) h'e[rM- 

(2) h"{s) + dist{h'{s),n) > i. 

(3) EiT < h"{s) < e.r + 1 for se[i,i + l]. 

(4) \h"'\ < h". 

The proof of the lemma is fairly straightforward and uses only the fact that 
{Ei) is slowly varying and summable. The second part is needed in order to 
avoid the eigenvalues of the Hermite operator. 

Lemma 6.3. Let t, (Eij) and (Ei) be as in Lemma \6. 1\ Then there exists a 
smooth spherically symmetric function 

: M X M" ^ M 

with the following properties: 

27 



(1) (Bounds) The function (j) is supported in \y\ < Ir^^'^ and satisfies 

(6.10) < <P{s,y) < e,T, \d,<P{s,y))\ + \d'Ms,y)\ < e^T 
and 

3 

(6.11) 5^ (1 + \y\)'\Dl+'di<j>\ < 6,ri/2 for t < s < z + 1, 

l,k=0 

(2) (Monotonicity) 

(6.12) dr(p{s,y)^ey^ for {s,y) e Aj, {i,j)eA{r), j>jii) + l 

(3) (Convexity) 

(6.13) il + \y\)d^,(P{s,y)^e,,T"^ m A„ (2,j) G ^(r). 
Proof. Let 

0,(y) = Ve2^- + \y\\ J > 0. 
We fix a smooth partition of unity 1 = "^rjls — i) and define 

In aj{s) = r](s — i) InSij. 

i 

These functions satisfy the bounds 

(6.14) aj{s)K,eij, i<s<i + l, la^-l, |aj|, |aj'| < a^. 
Their sum satisfies 

[lnT/2]+2 

a{s) := o-jis) ^ ^i, i < s < i + 1. 

j=0 

We define 

[lnr/2]+2 

<P{s,y) = Th{\y\r-'/') E ^MM\y\) 

j=0 

where x is a smooth function supported in [—2, 2] and identically 1 in [— |, |]. 
We verify the properties: 

< 0(s, y) < a(s)r < e^r, i< s <i + l. 

The remaining part of (16.101) follows from (16.141) . Estimate (16.111) is a conse- 
quence of (16.141) and 

|(l + |y|)^D'=+V,(y)l<l, 0<A;<3. 

The upper bound from (I6.12p is covered by (16. lip and the lower one follows 
from 

in Aij with j > j{i) where we use ~ Ei. The assertion (I6.13P follows from 
immediate bounds on second derivatives of the 0j. □ 

28 



Our aim in this section is to prove Li^ Carleman estimates for the variable 
coefficient operator P with the exponential weight 

e^^ 4 't^'e — sr. 

Here 

= h{s) + 62(p{s,y). 

where 82 is a small constant and h and are as in in Lemma 16.21 and 16.31 

The calculations are involved. For a first orientation we outline the key part 
of the argument for the constant coefficient heat equation. Using the change 
of coordinates of Section [3] we transform the problem to weighted estimates for 
the operator Pq = dg + H and the exponential weig ht e^^'^y\ This translates 
to obtaining bounds from below for the conjugated operator 

Lemma 6.4. Let r be large enough. Let ip he as in (16.91) with h, as in the 
above two Lemmas \6.^6.3\ with 62 ^ I ■ Then the operator Pq^^ satisfies the 
bound from below 

\\ih")h\\' + 52T-\\\al,\/v\\' + WalW^vf) < \\Po,^v\\l. 

for all functions v supported in < 9r} where the weights amt, (i± are 

defined by 

Proof. We decompose Pq^^ into its selfadjoint and its skewadjoint part 
where 

(6.15) LS,^ := -A, + y'-^s- V^,^ := 5, + + d^y. 
Expansion of the norm gives 

(6.16) \\{Ll^ + Ll^)v\\l. = \\U,^^v\\l, + \\Ll^v\\l. + {[U,^^,LlJv,v) 
The conclusion of the lemma follows from the commutator bound 

(6.17) {[L^o,^,Ll^]v,v) > \\{h'f^vr + 62T-\\\al,Vvr+\\alVM') 
The commutator is explicitly computed 

Since ^2 <^ 1 the ffist term has size h"{s). The second one is nonnegative since 
is convex for < 9r. 

The Hessian of the radial function can be written in the form 

(6.18) Ipyy = iPrrTl ® Tl ^ \ - Tl ® Tl ] 

\y\ \y\ r V \y\ \y\J 
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One can see that the radial and angular derivatives carry different weights. 
Our construction of (f) guarantees that 

hence the weight ilj^r can be used for all derivatives. For the size of the two 
weights we have 

/ ~ 4 > 4 

r 

This gives the last two terms in (I6.17p . 

It remains to see that the remaining terms in the commutator are negligible 
compared to the first term on the right hand side of fl6.17p . For this we use 
the bound (16.111) to conclude that in Aij we have 

I - Ay^y + Ai,yi,,y " A^^] < 5 26 iT < 62h" 

□ 

To switch to operators with variable coefficients it is convenient to extend 
the weights to the full space and to regularize them. Precisely we shall assume 
that 

^4 ^„ „A , _ A :c l„.|2 

(6.19) 



0'intis,y) ^ SiT in Aij if \y\ > r 

(s, y) ^ Eijil + in A^j if < r. 



4 

Observe that the two cases above match since Ei ~ Sij in the region where 
y"^ ^ T. We also introduce a modification a of amt which is used to include the 
effect of the spectral gap in regions where we have very little convexity: 

a'^{s, y) ^1 + EiT in Aij if \y\^ > t 

(6.20) ^ _^ 3 

a^{s, y)^l+ Eij{l + \y\) in Aij if \y\^ < r. 

Finally we choose a±_ with the properties 

supp a_L C \J{Aij : - I < j Inr} 

(6.21) a\_{s,y)<E,{l + \y\)-'Ti in A,, 

3 1 

A.{s,y) ^ei{l + |i/|)~V2 in Aij if j(i) < j < - Inr - 1 

The bounds for the weights from above are assumed to remain true after 
applying powers of the differential operators 9^, dy and ydy to them. 

Consider now a the more general class of operators P with real variable 
coefficients given by fll.7p . We repeat the change of coordinates and write in 
the (s, y) coordinates: 

4e-^^P = - 2y^ + d.g'^d, + yd'^d, + did'^y, + y.e'^y, 
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This further leads to 

2 2 

where P is given by 

We rewrite it in the generic form 

P = Pq — ddd — ydy — ydd — ddy 

with PQ = ds + H. 

To simphfy as much as possible the proof of the main Carleman estimate 
we introduce a stronger condition on the regularity of the coefficients: 

\d\ + {y){\dy\+T-'^\dyy\ + T-^\dyyy\ + T--2\ds\) <5ieij inAij 

(6.22) J ^ 

\d\ + {y)i\dy\+T—2\dyy\+T-^\dyyy\+T—^\d,\) < 5l 

This improved regularity will be gained later on by regularizing the coefficients. 
We are now in the position to formulate the Carleman estimate. 

Proposition 6.5. Let r be large enough and (5i <^ ^2 ^ 1. Let ip he as in 

(16. 9p with h, (p as in Lemmas \6.2^6.3[ Assume that the coefficients g — In, d 
and e satisfy (16.221) . Then the following L"^ Carleman estimate holds for all 
functions u supported in {y < 9r}.- 

(6.23) 4( r-'2\\a^e'^D^u\\^T-^^\\a\e'^D^u\\\ <\\e'^Pu\\. 

\i=o,i,2 / 

Proof. After conjugation 

we decompose P^ into its selfadjoint and its skewadjoint part 

which for y'^ < 9r can be expressed in the generic form (see also (I6.15P ): 

= LS,^ + ddd + Td 

L), = Ll^ + T^dd,+d,d) 
with d satisfying (I6.22p . Then (I6.23P follows from 

(6.24) Yl r-^iS2\\al,D^vf + {h"fw^vr) + 6,T-'\\alVM' < ll^^^f • 
i=o,i,2 

The proof will consist of three steps. 
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Step 1: First we show that for v supported in < 9r} we have 
(6.25) 

T 

We compute 

\\P^v\\l, = \\U^v\\l, + \\Llf,v\\l, + {[L:f,,L^]v,v) 

We expand the commutator 

[L;, Ly = [LS,^, V.J + [ddd + Td, LIJ + r5 [l;, dd, + d,d] 

The main contribution in fl6.25p comes from the first commutator, for which 
we use (16.171) to obtain the terms on the left side of (I6.25p . 
The second commutator is estimated by 

m'',LlJv,v)\<5,i\\alM\' + r-'\\al,Vv\^ 

and the second term on the right is neghgible since 6i <^ 62. Indeed, we write 

[ddd + Td,Ll^] = -dkq'''di + r 

where the coefficients q, r have the generic form 

q = ds+ IPydy + ijjyyd + dljjyy, T = OdOAlp + T{ds + IPydy) 

Using the bounds (I6.22p for d and (16. lip for (p we estimate 

Finally, the third commutator is estimated in a similar fashion. We write it 
in the form 

T^L'^,dd, + d,d] = -dkq'''di + r 
where the coefficients q, r have the generic form 

q = T^{dy + ddy), r = T^^Ady + ddddy + dipys + dipyyipy + rddy) 

Using (I6.22P and (16.111) we obtain the same bounds for q and r as in the 
previous case. This concludes the proof of (I6.25p . 

Step 2: We use an elliptic estimate to show that for functions v supported 
in ^ 9r} we have 

(6.26) S2(^^r-^\\al,D^vr + T-V±D^v\\'^^^^ < 11^' 

The elliptic bound 

P^^ll +^11^11 ^ ^'^\\Dv\\ + ||(-A - h'{s))v\\ 
can easily proven by a Fourier transform. It implies 

\\D^v\\ + r\\v\\ < T^\\Dv\\ + \\{H - h'{s))v\\ + \\y'^v\ 
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We can replace H — h\s) by due to the pointwise estimate 

|(L; -{H- h\s)))v\ < 5,{\D^v\ + T^Dv\ + t\v\) 

Then (16:261) follows from (16:251) . 

Step 3: Here we use the spectral gap condition to improve our bound when 
h" < 1 and show that dOe]) imphes (l6^ . It suffices to show that if h"{s) < \ 
then 

\\v\\l2 + t-^DM < W^^h^ 

Indeed, let s G [2,2 + 1] so that h"{s) < |. Then h' has a positive distance 
from the integers. Also ^ 1 which implies that at time s we must have 

\9 - In\ < 6,T-\ \Dg\ < 5,r-\ < 52T-1 

Hence we may think of as a small perturbation of H — h'{s) and compute 



WW + T-^\\D\\\ <\\iH- h'{s))v\\ < \\Lrv\\ + (5i + 6l)\\v\\ + 5it-^\\DW 



where the last two terms on the right are negligible compared to the left hand 
side. The proof of the proposition is concluded. 

□ 

We want to reformulate the previous result in a more symmetric fashion. 
To do this we weaken the estimates slightly by using a coarser partition of the 
space. We distinguish three cases for i corresponding to the value of in 
Lemma 16.11 (v) . 

Definition 6.6. We define the partition Bij as follows. 

(1) If = (which corresponds to si > r^2 ) we set 

Bij = Aij, a, b^^a± 

(2) IfO < jii) < [lnr/2 + 2] (which corresponds to r^^ ^ ^ J we set 

Bij = Aij, b^a, b± = a± j > 
respectively 

Bio = \^ Aij, b^ a|^^^.(^j b± = on Bio 

j<j{i) 

(3) If j{i) = [lnr/2 + 2] (which corresponds to r^^ ^i) we set 

[lnr/2]+2 

Bio = [J Aij, b = 1, b_L = on Bio, 

j=0 

Heuristically the definition of the Bij partition is motivated by the fact that 
in regions Aij with j < j{i) the weight is ineffective, i.e. it changes by at 
most 0(1). Thus the convexity there is useless, and instead we rely directly 
on localized bounds for the Hermite operator. 
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Since the Eij are slowly varying b < a and bj_ < a± and we may replace the 
a's by 6's in the above proposition. 

To provide some bounds on the size of b and b± we introduce a function 
1 < r{s) < T2 which is smooth and slowly varying on the unit scale in s so 
that 

r{s) ^ e^'(*) s e + 

This describes the region where b is tapered off and b± = 0. Precisely, consider 
two cases corresponding to the three cases above. 
(1) If r(s) ^ 1 then we have the bounds 

Mr(l + r)-l< b\r, s) < Mr(l + r)-' 

{^■27) 

fei(r,s) <Mr(l+r)-i 

where the parameter M > 1 is defined by M ^ e(s)r^. 
(23) If r(s) > 1 then 

28) rr(s)"5(r(s) +r)"i < b\r,s) < Tr{sy\r{s) + r)'^ 

with approximate equality when r < r(s) and approximate equality on the 
right when r = t2 . 

By slightly changing b and b± we may and do assume that the functions b 
and b± are smooth with controlled derivatives. Thus b and b± are smooth on 
the unit scale in s and on the dyadic scale in y, and their derivatives satisfy 
the bounds 

(6.29) \bs\ + {rs + r)\br\ + {rs + rf\brr\ <b, < 9r 
and 

(6.30) \b^ + {r, + r)\b^,\ + {r, + rf\b^,,\<b^ + b < 9t 
In addition we have 

(6.31) supp bj_r C {r > r.,} 

Using the functions b and b± we define the Banach space with norm 



\v\\h = \M\h + ^''^'\KDiv\\h 
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Then the symmetrized version of Proposition 16.51 has the form 

Proposition 6.7. Assume that the coefficients of P satisfy 06.221) . Let ip he 
as in (16.91) with h, as in Lemmas I6'.i^6'.^ Then the following Carleman 
estimate holds for all functions u supported in {y < 9r} .■ 

(6.32) l|e^(^'^)«||^o<||e^(^'^)p«||(^o). 
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Proof. Conjugating with respect to the exponential weight, the bound (16.321) 
is rewritten in the form 

(6.33) \\v\\xo<\\P^v\\^xor 
Observing that 

we introduce the operator 

Q = Q{\y\, (-As,.-0^), qir, A) = {b\r) + r-'T-X{r)X'f^ 
Then the inequahty (16.331) can be written as 

WQvWl^ < \\Q-'P^v\\l2 
whereas inequahty (16.231) imphes 

Hence it is natural to apply (I6.23P to the function w = Q~^v, which solves 

P^w = Q"'P^v + [Q, P^]w 

Thus (16.331) would follow provided that the commutator term is small, 

\\Q-'[Q,P^]w\\l2 < \\b^w\\L2 + r-'^^\\b^Vw\\L2 + T-'\\b^Dlw\\L2 

Unfortunately a direct computation shows that the smallness fails when j is 
close to even in the flat case, i.e. with P^ replaced by 

Po,v> = ds- A + y^ -^s-i^l 

To remedy this we introduce an additional small parameter 6 and use it to 
define a modification Qs of Q. We modify r(s) to rs{s) defined by 

r5(s)-2 = 58r(s)-2 + 5V-^ 

and use it to define the function 

bs{r,s)^ = 6~''T{r^ + rs{sy)-' 

We can still compare it with b, 

bs{r,s)'<6-%{r,sy 

Its usefulness lies in the fact that it is larger than b exactly in the region where 
the commutator term above is not small. 
The modification Qs of Q has symbol 

qs{r, s. A) = g(r, s. A) + bs{r, s) = {b^{r, s) + r"^r"^6^(r, s)A^)2 + bs{r, s) 
which satisfies 

(6.34) q<q5< S'^q 
We claim that it satisfies the bound 

(6.35) \\Qf[Qs,PM\L^^{S + c{6)6,) J2 r-'\\b'D^w\\L2 

j=0,l,2 
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Suppose this is true. Then we fix 6 sufficiently small, and for 6i small enough 
we apply fl6.23p to w = Qj^v. By fl6.35p he commutator term in the equation 
for w can be neglected, and we obtain 

which by f l6.34p implies that 

\\Qv\\l^<S-'\\Q-'P^v\\l2 

It remains to prove fl6.35p . 

I. We first calculate the commutator in the flat case, i.e. with replaced 
by Po,ip- Due to the spherical symmetry the only contribution comes from the 
radial part of the Laplacian and the s derivative. Hence using polar coordinates 
we compute 

Qs [Qs^ Po,ip] = Qs [Qsrr H ' Qsr + '^Qsrdr — QssJ 

Then is suffices to verify that on the symbol level we have 

(6.36) \qsrr\ + r-^\q5r\ + + r^^^lqsrl < 6qsb\l + r'V-^A^)! 

1.(1). We begin with the q component of qs. Using f l6.29p and fl6.30p one 
obtains 

\qrr\ + r~^\qr\ + \qs\ + T^/^|g.r| < {r{s) + r)~^T2q 
Thus it remains to show that 

(r(s) + r)~V^g < %6^(1 + r"V"^A^)5 

Optimizing with respect to A it suffices to consider the cases A = respectively 
A = rr2, where the above inequality becomes 

(r(s) + r)-V^(6 + b^) < 5{bs + b)b^ 

or equivalently 

{b + b^)bl<5{bs + by 

which is true since by (16.270 and (16.281) we have b^bi < b^ while bi < Sbs- 

1.(2). Next we consider the b^ component of g^, for which it suffices to prove 
that 

(6.37) \bsrr\ + r-^\bsr\ + \bss\ + r^'^\b8r\ < Sbsb^ 
I. (2). (a). For the s derivative we compute 

bs r'^ + rg{s) 

therefore we want to show that 

ir'sis))s<S{r' + rj{s))b' 

We optimize the right hand side with respect to r. The minimum is attained 
when = min{r^(s), r}. We need to consider two cases: 
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I.(2).(a).(i). If rs{s) < then r|(s) ^ 6^^r'^{s) and r > 5^^r'^{s). Hence 
using the estimate from below in (16.271) and (I6.28P we obtain b'^{rs) 
Then the above bound for r = rs follows since \{rs{s)~'^)s\ < rs^s)^"^. 

L(2).(a).(ii). If rs{s) > A then (ESHD we evaluate h'^{A) ^ T^r{s)-^. 
Then the above bound becomes 

Since < r(s)~^ it suffices to show that 

The worst case is r(s)^ = 6^t, rs{s) = when it is verified directly. 
I.(2).(b). For the r derivatives the last term is the worst. Since 

bs r2 + r|(s) 

we want to show that 

Optimizing with respect to r the worst case is when = min{r|(s), r}. 

I.(2).(b).(i). If r|(s) < r then rs{s) ~ 5~^r{s) therefore for r = rs{s) the 
above relation becomes 

r^<rV(s)62(r5(s)) 
which follows from the bound from below in fl6.27l) and fl6.28l) . 

I. (2).(b).(ii). If > r then as before we evaluate 6^(r5) ^ r3r(s)^5 and 
rewrite the above bound as 

r < 5r|(s)r4r(s)^2 

The right hand side is smallest either when rs{s) = r3 and r(s) = 5^r^ or 
when rs{s) = 6^^t2 and r{s) = ra. In both cases the inequality is easily 
verified. 

II. Now we deal with the general case, which we treat as a perturbation. 
Since we do not care about the dependence of the constants on 6 to keep the 
notations simple we include bs in b and work with Q instead of Qs- Thus in 
the computations below we allow the implicit constants to depend on 6. 

Suppose that A is a pseudo differential operator of order 1 and let r] be any 
Lipschitz function. Then 

(6.38) \\[A,r^]fU.<\\f\\L^. 
We write 

g(A) =b+{b^ + ^7r)3 - b =: b + gi(A). 
Even though qi has order i, we treat it as an operator of order 1 and estimate 

sup(A)'=-i|gi('=)(A)|<-^ 
A ror 2 

^the equality holds on the right when r'^ ~ t 
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Hence, for each r we obtain the bound on the sphere S"' ^ 

(6.39) ||[g,^]/||^,<il_||,^||^,^(5„_,)||/||^,. 
As a consequence, it also follows that 

(6.40) \\[Q,V^mL- < -^ll^ll W-)I|V^,/I|l- 

ror 2 

where Ve stands for the vector fields Xidj — xjdi generating the tangent space 
of ^"-1. 

To use these bounds we write the difference — in polar coordinates, 

P^ - Pi = P'.dl + P^dr + P^ 

where Pg are spherical differential operators of order j. Modulo zero homoge- 
neous coefficients which are polynomials in xr~^ we can write 

P^ = d, P^ = dr-^Ve + rU + dy 

p2 = dr~^Vl + {rh + dy)r~^Ve + rd + r^dy 

where d stands for coefficients satisfying (16.221) . In the support of b± we have 
a ^ b therefore our regularity assumptions on d show that for fixed r we have 

The coefficients involving dy satisfy better Lipschitz bounds and are neglected 
in the sequel. 

We expand the commutator 

[g, P^] = IQ' + Pe^Qrr + IQrdr) + PJQ. 

j=0,l,2 

Using the trivial b~^ bound for and (16.171) . (I6.4UI) we estimate the first 
term, 

\\Q-\Q.Pl\dl-^w\\^.<b-Hb\-^^_ J2 r'-HD^w^,. 

j=0,l,2 ^"''"^ i=0,l,2 

This is bounded by the right hand side in (I6.35P since 
The second term in the commutator is estimated by 

WQ'^PeiQrr + 2Qrdr)w\\L2 < ^ ( || Qrr^ || L2 + || Qr^r^ || L2 ) 

This is bounded by the right hand side in (I6.35P provided that 

\qrr\ + T-^qrl < 6^(1 + T^^r-^A) 

which follows from (I6.36p . The third term in the commutator is treated simi- 
larly. This concludes the proof of the proposition. □ 
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To conclude our study of the Carleman estimates we need to also pay 
some attention to elliptic estimates. The conjugated operator is elliptic in 
the region {y^ + > 4r}. Precisely, in this region we have the symbol bound 

\L;,{s,y,0\>y' + e 
Consequently, we can improve our estimates in this region. We consider a 
smooth symbol ae{y,C,) with the following properties 

supp tte C {y^ + > 8r} 

We define the space X2 with norm 

(6.41) \\vfx, = \\v\\lo + \\a:iy,D)vr 
The dual space X| has norm 

(6.42) 11/11^. = inf{||M|^^o), + II/2IP; / = /i + a:{y,D)f,} 

We note that due to the elliptic bound for high frequencies, we also have the 
dual bounds 

(6.43) T-'^\\bDv\\ < \\v\\x2, llV/IUi < lirVll 

Then our final Carleman estimate is 

Theorem 5. Assume that the coefficients of P satisfy ( I6.22p . Let ip be as 
in (16. 9p with h, as in Lemmas \6.2)i6.3[ Then the following Carleman 
estimate holds for all functions u for which the right hand side is finite: 

(6.44) \\e^^'^y\\\x:, < We'^'^'^y'^PuWx* 



Proof. We first prove the result using the stronger assumption (I6.22p on the 
coefficients. After conjugation we have to show that 

(6.45) \\v\\xo < \\P^v\\xi 

We consider two overlapping smooth cutoff symbols Xi = Xiiv"^ + ^^) and 
Xe = Xeiy'^ + i"^) ■ The interior one Xi is supported in {?/^ + ^^ < 7r} and equals 
1 in {y^ + > 6r}. The exterior one Xe is supported in {y^ + > 4r} and 
equals 1 in {y"^ + < 5r}. We need the following bounds for Xi and Xe'- 

Lemma 6.8. a) The operator Xi{x , D) satisfies the bound 

(6.46) \\x^ix,D)f\\^xoy<\\f\\xs 

b) The operators Xj(x, D) and x^{x,D) satisfy the following commutator 
estimates: 

(6.47) \\b-'[x^{x,D),P^]v\\<r--^\\bv\\ + \\xev\\ 



1 



(6.48) \\[Xe,PM^^~'\M 

39 



Proof, a) By duality the bound fl6.46p is equivalent to 

\\Xi{x,D)v\\x2 < ll^^llxo 

We have 

\\ae{x,D)x^{x,D)v\\<r-''\\f\\ 

since the supports of the symbols (1 — Xe{x,^)) and ae(x,^) are 0(r3) sepa- 
rated. Then it remains to show that 

\\Xi{x,D)v\\xo < \\v\\xo 

which is fairly straightforward and is left for the reader. 

b) We now consider the bound (16.471) . Commute first Xi with ds + H — h'{s). 
We have 

[Xi{x, D),ds + H- h'{s)] = [x^{x, D), H] 

Since Xe = 1 in the support of Vx,^Xi the Poisson bracket of Xi and + 
vanishes, by standard pdo calculus we obtain 

II [x.{x, D),ds + H- h'{s)]v\\ < \\xev\\ + r-^||^i;|| 

The difference — [dg + H — h'{s)) can be expressed in the form 

P^-{ds + H- h'{s)) = dgd + T"^{gd + dg) + rg 
where the function g satisfies the bounds 

1^1 + {y)\9y\ + {y)\gyy\ < £i 

These lead to an estimate for fixed s G [i, i + 1], 

\\h-\x^{x,D),P^ ~{ds + H- h\s))]v\\ < e,T^2\\{y)-^b-'v\\ 
Then (16.471) follows since 

Finally, the proof of the estimate (16.481) is similar but simpler. 

□ 

We continue with the proof of the proposition. For the nonelliptic part we 
apply (16.321) to the function Xi{x,D)v which is supported in {y"^ < 9r}. This 
gives 

\\Xt{x,D)v\\xo < \\xi{x,D)P^v\\(^xoy + \\[xi{x,D),P^]v\\L2 

For the first term on the right we use the bound (16.461) while for the second 
we use (16.471) . This yields 

(6.49) Wxiix, D))v\\xo < \\P^v\\x', + ||6t;|| + ||xet^|| 

On the other hand for the estimate in the elliptic region we compute 

(6.50) {{x:rv, p^v) = {x:v, + (x>, [x:, p;]v) 
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For the first term we split into H — h' plus a perturbation. Using pointwise 
bounds for the coefficients of we obtain 

\U^v -{H- h')v\ < 5i((r + y^)\v\ + t^Dv\ + \D\\) 

which shows that 

(X>, L;x» = {x:v, {H - h')x:v) + 0(5i(x>, {H + r)x») 

The symbol oi H — h' is elliptic in the support of Xe, therefore a standard 
elliptic argument yields 

(X>, {H + r)x» < (x>, {H - h')x:v) + Ct-'^M' 

for a large constant C. This further gives 

{x:v, {H + r)x» < L^xr^^) + Ct-'^M' 

Returning to (16.501) . we obtain 

We use (16.481) and then the Cauchy-Schwartz inequality to obtain 

(X>, {H + r)x» < II (i/ + ry'^P^vW' + r^t ||6t;f 

The first term on the right is properly controlled due to the straightforward 
estimate 

ll(^ + r)-^/|| <||/IU| 

Hence combining the above inequality with f l6.49p we obtain 

\\x\x,D))v\\x^^ + \\{H + Ty2x':{x,D)v\\ < \\P^v\\x* + r-^\bv\\ + \\x:{x,D)v\\ 

The last two terms on the right are negligible compared to the left hand side, 
therefore we obtain 

(6.51) ||t;|Uo + ||(if + r)^xr(a;,/^)^;|| < \\P^v\\x* 

Then (16.451) follows since Xe = 1 in the support of Og. 

It remains to show that the assumption (16.221) on the coefficients for (16.71) 
can be replaced by the weaker condition (11.151) . This is a direct consequence 
of (I6.43P combined with the following regularization result: 

Lemma 6.9. Let d be a function which satisfies (11.151) . Then there is an 
approximation gi of it satisfying (16.221) so that 

\9-9i\<b'r-' 

Proof. First we transfer (11.151) to the {s, y) coordinates. A short computation 
yields the equivalent form 

(6.52) IMIU-{A,,) + e^\\d\\Lipy(A^,) + lkllcr^(A,,) ^ 
where the new continuity modulus fhij is given by 

TrLij[p) = p + e 3 p3 
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Within Aij we regularize d in y on the 6y = t ^ scale and in s on the 
5s = e5r~3 scale, 

These localized regularizations are assembled together using a partition of unit 
corresponding to A^K In we compute 

\d-di\^ ^ij{(i~^6y + mij{5s)) ^ eij{e~h~^ + e"5r"^) < ef^.e"5r"^ < feV"^ 
while 

\dsgi\ < — ~ ^ij^ 

The bounds for higher order derivatives of gi follow trivially due to the fre- 
quency localization. □ 

□ 

7. W Carleman estimates for variable coefficient operators 

The variable coefficient counterpart of Proposition |4] uses the more convex 
weights constructed in Section [6l For convenience we write it in the (s, y) 
coordinates. Let r >> 1 and B{t) be as in fOHD . fOHD and ffL27j) . 

We define the function space X through its norm 

(7-1) \\v\\x ■= \\v\\x2 + ||^||/2(B(r);Lt°°L2) + || || i2(B(r);Lf'L|) 

where {p,q) is an arbitrary Strichartz pair, with X2 as defined in (16.411) . 
Its (pre) dual space has the norm 

(7.2) ll/IU* = , M , . Wf^Wxi + II/2IILIL2 + WMlp'li' 
^ ' j=ji+h+h 

Then we have the following improvement of Theorem [51 

Theorem 6. There exists ip as in (16.91) with h and (p in Lemma \6.2\ and 
l6'.o/i Then the following estimate holds for all compactly supported sufficiently 
regular functions u. 

(7.3) \\e^u\\x < \\e^Pu\\x* 



The relation between ip and the partition B{t) remains a bit mysterious at 
this level. If we replace it by the empty partition then the statement remains 
true for all ip with h and (p as in Lemma [6. 2[ The same is true for a partition 
into time slices of size 1. The convexity properties of (p allow a localization 
to the finer partition Bij as in 16.61 (and, as we shall soon see, to an even finer 
partition), q It is possible to choose (p and h so that the partition (Bij) is 
finer than the one defined by B{t). We assume in the sequel that ip has been 
chosen with these properties. 
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Proof. As usual this is equivalent to proving a bound from below for the con- 
jugated operator, 

(7.4) \\v\\x < \\P^v\\x* 

The main step in the proof is to produce a parametrix for P^. The key 
point is that the parametrix is allowed to have a fairly large error. This 
is because errors can be handled by Theorem O The advantage in having 
a large error is that it permits to localize the parametrix construction to 
relatively small sets, on which we can freeze the coefficients and eventually 
reduce the problem to the case of the Hermite operator. The properties of the 
parametrix are summarized in the following 

Proposition 7.1. a) Under the assumptions of the theorem there exists a 
parametrix T for P^ with the following properties: 

(7.5) l|T/||x<||/||x. 
and 

(7.6) \\P^Tf-f\U^<\\f\\x^ 
b) The same result holds with P^ replaced by P^. 

We first use the proposition to conclude the proof of the Theorem. Let 

P^v = f + g, ll/IUj + IMl^l^+lp'li' ~ WP^x* 

With T as in part (a) of the proposition we set 

w = v-Tg, P^w = f + g- P^Tg 

By (17. 5p we can bound Tg in X, therefore it suffices to bound w in X. On the 
other hand by (17.61) we obtain 

\\PMx', < + \\9- Pi.Tg\\x,, < \\P^v\\x* 

It remains to show that 

Ikllx < IIP^wllxj 

By Theorem [5] we can estimate the X2 norm of w and replace this with the 
weaker bound 

(7.7) \\w\\LooL2nLPLi < \\w\\x2 + ll-P^^lUa* 

This is proved using a duality argument and the parametrix T for P^ given 
by part (b) of the proposition. For / G X* we write 

{w,f)={w,p;Tf) + {w,f-p;Tf) 

= {P^w,Tf) + {wJ-P;Tf) 
Using both (17. 5p and (17.61) with P^ replaced by P^ we obtain 

\{wJ)\<{\\PM\xi + \Mx^)\\f\\x^ 

and (17.71) follows. This concludes the proof of Theorem [61 □ 
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It remains to prove the Proposition I7.1[ 



Proof of Proposition 7.1. The strategy for the proof is simple: On sufficiently 
small sets we can approximate the problem by one with constant coefficients 
and the properties of the parametrix follow from Section HI We use a partition 
of unity to construct a global parametrix from local ones. We obtain errors 
from 

(1) Commuting cutoff functions with the operator. Hence the partition 
has to be sufficiently coarse. 

(2) Approximating the variable coefficient operator by constant coefficient 
operators. Hence the partition has to be sufficiently fine. 

To elaborate on this we define the notion of a local parametrix: 

Definition 7.2 (Local parametrix). Given a convex set B we call T a (B-) 
local parametrix for P^ if for all f supported in B 

(7.8) l|T/||x<||/||x., 



x* 



(7.9) IIW-/|U.< 11/11 



and T f is supported in 2B. 

If T is a parametrix and rj is supported on 2B, rj = 1 on B then ipT is a 
local parametrix, but with constants depending on the commutator of P^ and 
rj. Vice verse, if (Bj) is a covering, {rjj) a subordinate partition of 1 and Tj 
are local parametrices then 

j 

is a global parametrix, because fl7.8p is obtained by summation, and 
provided 



and its adjoint 



This is obvious for the part and has to be checked for the other part. This 
strategy of constructing local parametrices leads, if it is possible, to estimates 
which are stronger than in Proposition 17.11 and Theorem El because we may 
replace the function space X by t^X{Bj) respectively. t^X*{Bj). 

In the first part of the proof we study the localization, and in the second 
part we provide the local parametrices. 
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7.1. Localization scales. Here we introduce a localization scale which is finer 
than the Bij partition of the space, and show that it suffices to construct the 
parametrix in each of these smaller sets. Precisely, the sets Bfj introduced 
below are the smallest sets to which one can localize the estimates for the 
operator P^. The choice of their size is not yet apparent at this point, but 
will become clear in the very last step of the proof, where we estimate the 
commutator of with cutoff functions on such sets. We consider three cases 
depending on the size of €{. 

(1) If Ei < then we use Bio as it is. 

(2) If r^^ < < T~2 then we partition the set -8,0 into time slices B^q of 
thickness 

5s = 

(3) If r^^ < Si and j 7^ then we partition B-ij into subsets B^^ which have 
the time scale, radial scale and angular scale given by 

5s = br/, 5y = rh-/, 5y^ = rh;.^^ 

This gives a decomposition of the space 

R X R" = y Bfj 

We also consider a subordinated partition of unity 

Suppose that in each set Bfj we have a parametrix T^j satisfying (17. 5p and 
(17. 6p . Then we deffne the global parametrix T by 

We have by an iterated application of Minkowski's inequality 

\\Xijf\\p{L^L2+Lp'Li') ^ WfWl^L^+Lp'Li' 

and the dual bound 

II ''^Tl^jXijf\\L°°L'2nLPL'i ^ \\Tl'jXijf\\p{L°°L^nLPLi)- 

Hence (17.51) for T would follow if we proved that 

iij^^iu. < ii4ii'^^2 

where vfj = T^jXijf are supported in B^j. This is trivial by orthogonality for 
the component of the X2 norm. It is also straightforward for the elliptic 
part since the kernel of the operator in the following sense: Let Xo ^ C°°(]R) 
be supported in [—3/2,3/2], identically 1 in [—1,1]. We deffne Xe{x,^) = 
Xo{5-{x^ + e)/T^)- Then 

ii5Zx>j;-iix, < iix>jiipx. 
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and the adjoint estimate holds since the kernel of is rapidly decreasing 
beyond the r~2 scale, which is much smaller than the smallest possible spatial 
size for namely 5y > r~4. 

It remains to consider the angular part of the X2 norm, which is best de- 
scribed using the spherical multiplier Q appearing in Proposition 16. 7[ The 
symbol of Q is smooth with respect to A on the r^rb'^b^^ scale therefore its 
kernel is rapidly decaying on the angular scale 66 = T^^r^^b^'^b']_ which cor- 
responds to 6y± = T'2b\b^'^. But by ([3]) and (16.51) this can be no larger than 
r~8 which is again much smaller than the smallest possible spatial size for B^y 
Thus orthogonality arguments still apply. 

Then we have 

[Po, = 2x>y + Xti^ + x> 
We claim that the right hand side is negligible in the estimate. For this it 
suffices to verify that 

Ixj/'l < < bl, \x'i\ < bl 

The last relation is trivial. For the first two we consider three cases. 

(1) If j = and Ej < Ct^^ then we need no spatial truncation. We are 
allowed to truncate at \y\ > Cr^ to separate the elliptic region, though. 

(2) If j = and ej > Cr'^ then 

lAj/ \ ^ y \Ayy\ rvj "-' 

while 

(3) Otherwise, 



while 



\-yij\ < p-i \-yij I < p-2j 

\A.y \ r-~j ^ 5 \A.yy\ r-~j ^ 



4 = aJ, = e.,rie-^>Ce-^^»r 



The results can be summarized by saying that it suffices to construct local 



parametrices in the sets 



7.2. Freezing coefRcients. Our first observation is that restricting the result 
in the proposition to a single region B^- allows us to freeze the weights 6, b^ 
in the X2 norms. 

Next we are interested in freezing the coefficients of P. We consider the 
same three cases as above: 

The case Si < r^^ . In this case we are localized to 

Bio = [i,i + l]x 

and we have 

bio^l, Eij < 
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By f l6.2p the second relation leads to 

\d\ < 

Then using also (16.43^ we can estimate the terms involving d in the expression 
(E22]) for P^, 

(7.10) \\dddv\\x*+r\\dv\\x*+T'^\\{dd + dd)v\\x^ < \\v\\x. 

Hence without any restriction of generality we can assume that d = in P^, 
which corresponds to taking g = !„■ 

We also observe that in this case we have 

101 <1, 10.1 <r~^ 
Then we can also drop the component of ip. Finally, since 

\hss\ < 1 

we can replace h by its linearization at some point in the corresponding s 
region. 

Conclusion: It suffices to prove the result when d = 0, ip{y, s) = ts. 
We note that the separation of r from integers is no longer needed due to 
the localization to unit s intervals. 

The case < Si. In this case we are localized to a region of the form 
and we have 

The second relation leads to 

\d\ < e-^'«r-^ 

Then (17.101) is still valid, so we can assume again that d = in P^. 
We also observe that in this case we have 

101 < 1, 10.1 < e~'^'^ 
Then we can also drop the component of ip. Finally, since 



we can replace h by its linearization at some point in the corresponding s 
region. 

Conclusion: It suffices to prove the result when d = 0, ipiVj = ts. 
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The case r ^ < Ei. In this case we are localized to a region of the form 

B-^ = [so, So + r-3e5£^r2] x B{yo, r^^e^e^r^)^ l^^l ~ g^' 
and we have 

;2 ^ 3 _3 

Using (16.21) it follows that 

1 _i i 

\9{s,y) - g{so,yo)\ <elf ^e^ 

Arguing as before, this allows us to freeze d within B^^. However, we note that 
we are no longer allowed to replace c? by 0. 

Next we turn our attention to the weight function ip. First we have 

which allows us to replace h by its linearization in s at Sq. 

Secondly, we claim that we can replace by its linearization at yo- In the 
radial direction we have weaker localization but a stronger bound 

In the transversal direction we have better localization but a weaker bound, 

\(i>yy\ < SiAe'^. 

The first bound allows us to obtain the relation 

- 0'(yo,so)| < £,7rte-2 ^ 6^ 
Using also the second bound we can write 

{(pyiy, s) - (f)y{yo, So))dy = Vrdr + Z^±<9_L 

where the coefficients and are smooth on the B^ - scale and satisfy the 
bounds 

\^r\<r-hl, \yr\<r-hl^^ 
Conclusion: It suffices to prove the result when d = g{so,yo), 'ipiu^s) = 
Ts + cy, \c\ < eiT2. 

Additional simplification in the highly localized case. Given the above simpli- 
fications we need to work with a constant coefficient operator which has 
the form 

= -dt + H - T + ddd + cd, \d\<eij, \c\ < Si 

We diagonalize the second order part with a linear change of variables to obtain 

P^ = -dt + H~T + cd + Oieij)y'' 

We can freeze the last term at yo and add it into r. To deal with c we make 
the change of variable 

y ^y- {s- so)c 
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Then our operator becomes 

P^ = -dt + A-{y- c{s - so)f + r 

and the s — Sq terms are neghgible due to the s locahzation. 

Conclusion: We can assume without any restriction in generality that g = In 
and ip = Ts. 

7.3. The localized parametrix. We begin with the global parametrix K 
constructed in Section [51 Then we define the parametrix Tb in B by 

Tb = X2bK, B = 4 

and show that it satisfies (17.81) and (17. 9p . 

The U' part of (17. 8p follows directly from (15. 4p . It remains to prove the X2 
part, 

The elliptic part of the X2 bound, namely 

is obtained by an argument which is similar to the one beginning with (I6.50p . 
For the rest we consider two cases. 

i) If J = then _B is a ball, and we can use (15.41) directly with R = d. 

ii) If j > then B is contained in a sector B C -Bn.d but may be shorter 
than R. This is why we can use (15.41) for the angular part of the X2 norm, but 
not for the part. However, the part can be always obtained by taking 
advantage of the time localization, 

W'TBfU^ < WTBfh^L^ < \\fhrL^^L.'L.' 

It remains to consider the error estimate (17.91) . We have 

f-{d,-H + r)TBf = [X2B, ds-H + T]Kf = [x2b, - H + T]x,BKf 

But arguing as above x^bKoI satisfies the same X2 bound as X2BKof ■ Hence 
it suffices to show that 

[X2B.ds~H + r]:X2^X* 

This is where the dimensions of the set B are essential; they are chosen to be 
minimal so that the above property holds. We have 

[X2B, ds- H + t] = -dsX2B + {.dyX2B)dy + dy{dyX2B) 

= -dsX2B + idrX2Br)dr + dr{drX2Br) + {,d^X2By)d^ + 9_l(9j_X2b) 

For the first factor we use the bound 

\dsX2B\<b% 

For the radial derivatives of X2B we combine (16.431) with 

\drX2B\ < r'hl 
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Finally, for the angular derivatives we use the angular H 2 norm in X2 and the 
bound 

\drX2B\ < 

□ 

8. The gradient term 

In this section we consider the full problem, i.e. involving also the gradient 
potential W. Ideally one might want to have a stronger version of Theorem [6] 
which includes additional bounds for the gradient, more precisely for 

But such bounds cannot hold, for this would imply that one can improve the 
indices in a restriction type theorem. To overcome this difficulty we proceed 
as in |T7j, using Wolff's osculation Lemma. Wolff's idea is that by varying the 
weight one can ensure concentration in a sufficiently small set, in which the 
gradient potential term is only as strong as the potential term. Thus we still 
obtain a one parameter family of Carleman estimates, but with the weight 
depending not only on the parameter but also on the function we apply the 
estimate to. 

Given a gradient potential W satisfying (11.191) . we first readjust the param- 
eters Eij, Si constructed in Lemma 16.11 in order to insure that we have the 
additional condition 

||m^||l"+2(a[) < 

Then we begin with the spherically symmetric weights ip constructed in Sec- 
tion [6] and modify them as follows: 

(8.1) ■^{s,y)=i^{s,y) + 5k{s,y) 

where the perturbation k is supported in {|?/| < 9r} and is subject to the 
following conditions: 

(8.2) \d:dPyd{k{s,y)\ < sy-^ . e [^,^ + 1] 

Here 5 is a sufficiently small parameter. 

In order to prove the strong unique continuation result in the presence of 
the gradient potential W we need the following modification of Theorem ([6]): 

Theorem 7. Assume that (11.151) holds. Then for each r > and W subject 
to 

I|w^IIl"+2(A[) < 

and each function u vanishing of infinite order at (0, 0) and 00 there exists a 
perturbation k as in (18. 2p so that 

(8.3) ||e*M||x + ||e*iyVM||x* + ||e*V(W^u)||x-+r^||e*iyM|U- < ||e*(^-)Pn||x* 
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Here and in the sequel we will omit indices for W. After returning to the 
[x, t) coordinates and taking fll.l9p into account this implies Theorem [31 The 
reader should note that the choice of depends on both u and W. This is 
essential since for fixed (18. 3p cannot hold uniformly for all u and W. 

Proof. Up to a point the proof follows the steps which were discussed in detail 
before. We outline the main steps: 

STEP 1: Show that the Carleman estimate fl6.23p holds with replaced 
by \1' for all perturbations k as in (18.20 . The new conjugated operator is 
obtained from after conjugating with respect to the weig ht e'^^y^'l This 
adds a few extra components to the selfadjoint and skewadjoint parts, 

Ll=U^ + kl + ks + 2ky{ll)y + d) 
L\, =Llf,- ky{l + d)d - dky{l + d) 

Observing that we can write 

k^ + ks + 2ky{tljy + d) = Td, ky{l + d)=T^d 

with d as in (16.220 we conclude that the conjugated operator P^ retains the 
same form as P^, therefore the proof of (I6.23P rests unchanged. 

STEP 2: Show that the symmetric Carleman estimate (16.230 holds with 
ip replaced by for all perturbations k as in (18. 2p . Since Pq, has the same 
form as P^, this argument is identical. 

STEP 3: Show that the symmetric mixed fl Carleman estimate in 
Theorem [6] holds with ip replaced by for all perturbations k as in (18. 2p . 
Since Pij, has the same form as P^, this argument is also identical. 

STEP 4: Decompose W into a low and a high Hermite- frequency part, 

W = Wio^ + Whigh, Wio^ = x\{x, D)W 

where the smooth symbol x\{^iC) is supported in {x^ + < 81r} and equals 
1 in the region {x^ + < 64r}. Then we show that the high frequency part 
of W satisfies the desired estimates for all perturbations k as in (18.20 . namely 

(8.4) \\e''WhighVv\\x* + \\VWhighv\\x*+T^\e''Whighv\\x* < \\e''u\\x\\W\\Ln+2 
After conjugation this becomes 

\\Whigh^v\\x' + \\VWhighv\\x' + r^\Whighv\\x* < ||f IUI|W^||l"+2 

We only consider the first term on the left. The second one is equivalent 
by duality, and the third one is similar but simpler. We divide v into two 
components, 

where the smooth symbol xli^^O is supported in {x^ + > 9r} and equals 
1 in the region {x^ + > lOr}. 
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For the high frequency component of v we use the part of the X"^ norm 
to estimate 

WWhinhNYlvW 2(n+2) < ||W/iig/i||L"+2||VXe^^||L2 < ||iy||in+2||t;||x 
L n 

For the low frequency component of v it is still possible to estimate directly 
the high frequency of the output, 

\\xliW^,,,,\/{l~xl)v)\\H-^ <r-^\m,^Vil-xl)v\\L^ 

<\\Whigh\\L" + 2r~^\V{l - yI)v\\^ 2(n + 2) 
<||iy||Ln + 2||f II 2(n+2) 

L TL 

Finally, the last remaining part has a much better estimate, 

11(1 - xl){Wn,aHV{l - Xl)v)\\ < r-^||l^|U.+2 ||t;|| 

which is due to the unbalanced frequency localizations of the two factors. 

Due to the estimate (18.41) . it suffices to prove (18. 3p with W replaced by Wiow 
This allows us to replace the term V{Wiowv) by 

where we can estimate 

||VW^;o^|Un+2 <r5||iy||i„ + 2 



Hence without any restriction in generality we can drop the third term in (18.31) 
and show that we can choose the perturbation k so that 

(8.5) We'^WVullx* + T^\e''Wu\\x* < ||e*(")PM|U* 

STEP 5: Show that, given u and W, we can choose the perturbation k 
so that (18.51) holds. At this stage we no longer need the full X* norm for 

2(n + 2) 

the W terms, it suffices instead to consider the L "+4 norm. Begin with the 
unperturbed integral 

, , , 2{n + 2) , 1 , 2(n + 2) 

F^dxdt, = \e^WVu\^+^ + \T^e'f'Wu\^+^ 

We can select a subset / of R consisting of time intervals of length 1 with unit 
separation at least 8 so that 

n 

F,i,dxdt S I I F,i,dxdt 





By a small abuse of notation we label 

J = |J/„ hc[t~l,t + l] 
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We define a family of perturbations k depending on parameters bi, ai by 




Due to the choice of the intervals Jj it is easy to see that after changing the 
weight to \E' we retain the concentration to a dilate of /, 

F^dxdt < / / Fi^dxdt 
J 31 Jr" 

The choice of the parameters 6j, CTj can be made independently for each i. 
We consider two cases. 

i) Suppose Ei < r~2. Then the choice of the parameters is irrelevant since 
in 3/j we can estimate 

IIc'^W^VmII 2in+2) + T^\e^Wu\\ 2{„+2) < ||W^||Ln+2(||e*VM||L2 +r^||e*M||i2) 

<(r-V2||e'fv«||^. + ||e*n|U.) 

ii) Suppose ^ r~5. Then we need to choose the parameters 6j, cij in a 
favorable manner. This choice is made using Wolff's Lemma: 

Lemma 8.1 (Wolff's Lemma [29] )• Let fi be a measure in and B a convex 
set. Then one can find hk & B and disjoint convex sets C M" so that the 
measures e^*'=/i are concentrated in Ek, 



[ e^'^dfi >]-j e-'^d^ 

J El. J El. 



'I 



„+l Tl±l 



and 

We apply the lemma for the measures 

In our case we have 

Si = fe. ([-r,r] X B(0,r5)) , \B. 

Hence we can find parameters h\ and a\ and convex sets E\ C 3/, x 5(0, 3r2 
so that the corresponding measures F^^. are concentrated in E\ with 

V \E^\~^ > 5"+^^ 

At the same time we have 

~ \W\''^^dxdt<€l^^ 
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Hence we can choose k so that 



\W\''+^dxdt<e,T-'^\E^\-' 



which by Holder's inequahty leads to 



1 1 



(8-6) IIW^II . <£i^'r-^ 

Denoting this index k by kii) we can write 

||e*iy(V,r5)u|| 2(„+2) < ||e'^W^(V,r5)u|| 2(„+2) 

Decomposing the function v = e^u into low and high frequencies we further 
estimate 

||e*iy(V,r^)w|| ^j^<\\W{vM){l-x]{x,D)){e'u)\\^ 

PL 1+4 "+4 (li) 

The first term on the right is estimated as in Step 4, 
\\W{vM){l-x]{x,D)){e'u)\\ ^ 2(^ 

tfL "+4 [U) 

< \\W\\iooLu+2\\e^u\\x 
It is only for the second term on the right that we need to use fl8.6p : 

^ II * II 

5- e U\\ 2(n + 2) 

^ II * II 
^ l|e u\\x 

The proof of the Theorem is concluded. □ 

Appendix A. The change of coordinates 

Suppose that the coefficients g satisfy fll.151) . In this section we verify that 
we can change coordinates so that fll.lSp and (11.161) are both satisfied. Due 
to the anisotropic character of the equation we must leave the time variable 
unchanged and consider changes of coordinates which have the form 

{t,x)^{s,y), s = t,y = x{t,x). 
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The expression for the operator P in the new coordinates is 

P = dt + dkg'\t,y)di + ^d''di 

where the new coefficients g, d are computed using the chain rule, 

There is a price to pay for this, namely in the new coordinates we obtain 
lower order terms which cannot be treated perturbatively. Instead we obtain 
coefficients which have the same regularity and size as gi — In- 

The Lipschitz condition fll.151) ensures that g has a limit at (0, 0) so we 
assume that g is continuous. After a linear change of coordinates we may and 
do choose g with (7(0, 0) = J„. Again by (11.151) this implies 

(A.l) \g{t,x)-In\<t:l. 

Proposition A.l. Let g he a metric which satisfies (11.151) with (7(0,0) = In- 
Then there is change of coordinates {t,y) = {t,x{t,x)) which is close to the 
identity 

(A.2) < 1 

and has regularity 
(A.3) 

sup II (t + x'r'/\tdtnit + x'Y/'d,fx\\iHAiry,L^) « 1, 2<2a+ \f3\ < 4 

SO that in the new coordinates both functions g and d satisfy (11.151) . while g — In 
and d satisfy (I1.16P . 

Proof. Consider the covering of the [0, 2] x B{0, 2) = UAij with an associated 
smooth partition of unity rjij. We can assume that the functions r]ij satisfy 

(A.4) \d^d^,7^,,\<Co.pt~-{t + x'r^ 



We choose the points 



(ti,Xij) — (e ,e ~^'') G 



and insure that rjij = 1 near (t j , Xij ) . By (11.151) we have 

(A.5) sup ^ \g{ti,Xij)-g{ti,Xi^^j^i))\ + \g{ti,Xij)-g{ti+i,x^i+i)j)\ <^l. 
(i,j)eA{T) 

Within a fixed set Aij we consider the linear map defined by the matrix 

Xij 9 ^ {tiy -^ij ) • 

It transforms the coefficients at {ti,Xij) to the identity and has the desired 
properties within Aij. We assemble the maps defined by Xij using the partition 
of unity, 

X(tj x) ^ ^ Viji^^ X)XijX- 
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Then 

Vx(i,a^) - In = ^iVriij)xijX + r]ij{xij - /„), 
Let {t, x) G Ai^j^. Since ^ ^Vij = we have 

^xXit, X)-In= ^^Vijit, x){Xij - XioJo) + VijiXij - In)- 

riij{t,x)>0 

The first term on the right hand is small by flA.SP (for Xij) and the second one 
by flA.ll) therefore the smallness of Vx — In follows. 
For the second order spatial derivatives we write 

Dixit, ^) = Y D'^Vijit, x)xijX + 2D^T]ij{t, x)xij 

= Yl ^xVijit, x){Xij - XioJo)^ + '2D^r]ij{xij - XioJo)- 
Hence by flA.4p and (lA.Sp (again for Xij) "we obtain 

\\{\X\^ + tY/^Dlx\\lHA(r);L^) < 1- 

Also 

dtX{t,x) = ^5t?7ij(xij - Xiojo)^ 

gives the desired bound for the time derivative. A similar computation yields 
the bound for the higher order derivatives in ( 1A.3I) . 

Consider now the new metric g. Since both Dx and {Dx)~^ are Lipschitz on 
the dyadic scale with 1^{A{t)) summability, from fll.lSp for g we easily obtain 
fll.lSp for g. In addition, our construction insures that g(ti,Xij) = In- This in 
turn leads to the bound 

which shows that (11.151) for g implies (I1.16P for g — 1^- 

It remains to consider the lower order terms. From dtX we obtain coefficients 
d of the form 

d = tY ^tVijXij 

Within AjgjQ this gives 

d = tY 9tVij{t, x){xij - Xiojo) 

The functions tdtrjij (t, x) are bounded and smooth on the dyadic scale, while 
the 1^{A{t)) summability comes from the Xij ~ Xiojo factor due to (lA.Sp . Hence 
both (11.151) and (I1.16P are satisfied. 

The contribution of d'^x to d has the form 

d= %{dxX)~\dlx)9 

There is no singularity at x = since x is linear in x for x^ <^ t. Then from 
(IA.2P and flA.3l) we obtain 
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with added 1^{A{t)) summability inherited from d^x- This is better than 
f ll.l6p . and in effect this term can be included in W and treated perturbatively. 
The bound (11 .16^ is also easy to obtain from the similar bounds for g and 
derivatives of x- 

□ 
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